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II. 1 SUMMARY 

This section of the Final Report documents the Task II #1 "Cabin Atmosphere/ 
Cabin Wall Thermal Interface" unit study area investigation. 

The thermal scale modeling criteria applicable to radiation-conduction-convec- 
tion systems are derived. 

Detailed consideration is given to four possible scale modeling techniques. 
Modified Material Preservation 
Temperature Preservation 
Scaling Compromises 
Nusselt Number Preservation 

The most promising techniques (scaling compromises and Nusselt number preser- 
vation) were chosen for an experimental investigation. 

Full scale and 1/4 scale models of a configuration, in which radiation, con- 
duction and convection heat transfer effects are important, were built and 
tested. 

The full scale and 1/4 scale model temperatures were correlated for scaling 
compromises of mass flux and heat transfer coefficient preservation. 

A thermal math model was developed for use in conjunction with the Nusselt 
number preservation technique. A Nusselt number correlation was developed 
for the free convection tests using the 1/4 scale model data. However the 
investigation of the Nusselt number preservation technique was not completed 
due to a lack of time and funds . 

If was concluded that either preservation of mass flux or heat transfer coeffi- 
cient may result in adequate thermal similitude depending on the system being 
modeled. For manned spacecraft heat transfer coefficient preservation should 
give good thermal similitude for both steady state and transient conditions. 

The Nusselt number preservation technique is workable however it is costly and 
difficult to implement. 
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II. 2 INTRODUCTION 


Most spacecraft thermal scale modeling studies have dealt with systems involving 
only the radiation and conduction modes of heat transfer. These studies range 
from the verification of scale modeling techniques using relatively simple con- 
figurations (see for example references 1-4) to the thermal scale modeling of 
an actual spacecraft concept (reference 5) . Thermal scale modeling of manned 
spacecraft involves the convection as well as the radiation and conduction modes 
of heat transfer. The similitude criteria and discussions of possible scaling 
techniques have been presented in the literature (see for example reference 6) . 
However, no experimental verification of the scaling techniques has been reported 
for manned spacecraft. A thermal scale modeling of free convection in heated 
enclosures has been reported (reference 7). An 8 x 8 x 8 foot prototype room with 
a convector heater and a 1/4 scale model were built and tested. Radiation- conduct ion 
scaling criteria were used for temperature preservation in the 1/4 scale model, 
however, no attempt was made to modify the convective heat transfer and air at 
atmospheric pressure was used in both model and prototype. 

In order for thermal scale modeling to be a practical tool for manned space- 
craft applications effective thermal scale modeling techniques must be developed 
for the cabin atmosphere/spacecraft cabin wall thermal interface. 

This report describes the development and demonstration of practical thermal 
scale modeling techniques applicable to radiation - conduction - convection 
systems with particular emphasis on the cabin atmosphere/ cabin wall thermal 
interface. 
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II. 3 SCALE MODELING CRITERIA 

The thermal scale modeling criteria may be developed from the equations which 
govern the system thermal energy balance. The governing equations when written 
in nondimens ional form result in dimensionless groups of parameters. These 
dimensionless groups must remain invariant for thermal similitude to exist 
between similar systems. 

The heat conduction within the solid elements of the system is governed by the 
Fourier equation 

? C ' sr ^ ^ Z T fa ( 1 ) 

and the heat transfer boundary condition at the solid/fluid interface is given 
by 

h ~ 4(e,-VT)„= f<r(-r+- - 4 (*» ■?%)„ ( 2 ) 

If the heat conduction in the solid is two dimensional (in the plane of the 
surface) then equation (1) is incorporated in equation (2) which becomes 

-4 & ■ vr f l (3. 

These equations can be written in nondimensional form as: 

Heat conduction 



Boundary condition at solid/fluid interface 



or for two dimensional conduction 
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Thermal similtude within the solid elements of geometrically similar systems 
having identical radiative surface properties is achieved by keeping the 
following dimensionless groups invariant: 


o 


Heat Conduction: 



for two dimensional conduction 


o 


o 


Volume heat sources: 
Surface heat flux: 



Heat transfer to fluid: 



[Aa 


The invariance of this last dimensionless group ^ assumes that thermal 

similitude exists for the fluid elements of the systems as well as for the 
solid elements . 


The fluid elements of the system are governed by the energy, momentum and 
continuity relationships. The energy equation may be written (without 
internal heat sources , viscous dissipation or radiation) as 

Pr (AS * rr f j = 4 t V*T t <7) 

and the momentum (Navier-Stokes) equation may be written for incompressible 
flow with buoyancy effects as 

+&-Vr)--VP + j)V‘ir- P f ?e(T}-T})e 9 ( 8 ) 

The continuity condition for incompressible flow is 


V'tr^O 


( 9 ) 


These equations may be written in non-dimensional form as 


Energy 

(&) (f\ in + RePr ** ' ‘ *** 


(10) 
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Momentum 


(£)(jr) f +Re%%-Gr0 f e t 


(ID 


Continuity 


££• £& = O 


(12) 


Thermal (and dynamic) similitude within the fluid elements of similar systems 
then requires the following dimensionless groups to remain invariant: 


Reynolds number: Re - 


Prandtl number: Pr - 


Grashof number: Gr = 


Fluid/solid transients 


A &L. 

M 

f- 


A* 

: /£]/& 
\?clK 4 U 


The last dimensionless group when kept invariant preserves the 

relationship between the transient response of the fluid and that of the solid. 
However, for thermal scale modeling applications involving the manned space- 
craft cabin atmosphere-cabin wall thermal interface, the fluid (gas) transients 
have negligible effect on the spacecraft thermal response. Consequently the 
invariance of the mm, group is not required for thermal similitude 

in this case. 


The thermal scale modeling criteria for the cabin atmosphere - cabin wall thermal 
interface are summarized in Table II-l. 
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TABLE II-l THERMAL SCALE MODELING 
CABIN WALL THERMAL 

CRITERIA FOR CABIN ATMOSPHERE 
INTERFACE 

Cabin Wall 


Heat Conduction 

(k/LT 3 )* = 1 
o 

two dimensional conduction 

(k6/L 2 T 3 )* = 1 
o 

Volume Heat Sources 

^ y* - 1 

Surface Heat Flux 

q * = 1 
s 

Transients 

(kt/pCL 2 )* = 1 

Cabin Atmosphere 


Heat Conduction 

(k f /LT 3 )*= 1 

Reynolds number 

Re* = 1 

Prandtl number 

P * = 1 ' 
r 

Grashof number 

Gr* - 1 

^Transients 

(k f pc/k p f c f )* = 1 


*The cabin atmosphere transient effects on the wall are generally small 
and adequate similitude may be achieved without meeting this criteria. 
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II. 4 SCALE MODELING TECHNIQUES 

Four scaling techniques were considered for the cabin atmosphere - cabin wall 
thermal interface. These are the Modified Material Preservation, Temperature 
Preservation, Scaling Compromises and Nusselt Number Preservation Techniques. 
Table II-2 gives a comparison of these techniques. The Modified Material and 
Temperature Preservation techniques are straightforward attempts to meet the 
scaling criteria given in Table II-l. The Scaling Compromises technique attempts 
to achieve thermal similtude without completely satisfying the scaling criteria. 
The Nusselt Number Preservation technique uses the thermal scale model to 
experimentally determine the Nusselt number as a function of Reynolds and Prandtl 
numbers and uses these results in conjunction with a thermal math model to deter- 
mine the prototype performance. 

II. 4.1 Modified Material Preservation 

The material preservation technique used for radiation-conduction systems keeps 
the same materials in model and prototype and meets the scaling criteria by 
increasing the temperature and heat flux in the model. The scaling criteria for 
systems involving convective heat transfer may be met by a '-’modified material 
preservation technique" which accounts for the variation of gas thermal con- 
ductivity with temperature. This technique requires the scale model temperature 
and heat fluxes be increased to a greater extent than does the normal material 
preservation technique. The required model to prototype temperature ratio is 
given by 



The temperature requirements as a function of scale ratio are compared in Figure 
II-l for these techniques with air at 70°F in the prototype. The modified 
material preservation technique also requires the wall material to be changed 
(or in the case of two dimensional wall conduction, the proper selection of 
wall thickness) to meet the scaling criteria at the higher temperature. 
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(4) The use of "Nusselt Number Preservation 11 requires a verified math model for the radiation-conduction aspects of 
the system. 
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Since the Prandtl number does not vary substantially with temperature (or from 
gas to gas) and since the Reynolds number may be adjusted to the desired value, 
their preservation in the scale model is easily accomplished * For equal Grashof 
numbers in the model and prototype the scale model pressure must be increased 
such that 

(tp 3 (#*)' * (l*) * ci4) 

The scale model pressure requirements for simulating air at 70°F in the proto- 
type is shown in Figure II-2 as a function of scale ratio. 

The use of non metals (seals, insulation, etc.) in the scale model effectively 
limits the use of the Modified Material Preservation technique to scale ratios 
greater than about 0.5* 

II. 4. 2 Temperature Preservation 

The temperature preservation technique maintains equal temperature in the model 
and prototype and meets the scaling criteria by requiring a reduced thermal 
conductivity for the model materials. This technique has been used successfully 
for radiation-conduction systems at various scale ratios, and has considerable 
flexibiltiy for two dimensional conduction systems where the model material 
thickness need not be scaled geometrically. The application of this technique 
to systems involving convective heat transfer depends on the availability of 
gases with substantially lower thermal conductivities than that of the prototype 
gas. Figure II-3 shows the thermal conductivities of various gases as a function 
of temperature. The scale ratios which are possible, using selected gases to 
simulate air at 70°F, are depicted in Figure II-4. 


As with the modified material preservation technique, the preservation of Reynolds 
and Prandtl numbers in the scale model is easily accomplished. The preservation 
of the Grashof number in the model requires the model pressure to be set such that 



(15) 
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Figure H-2: PRESSURE REQUIREMENTS FOR PRESERVA TION OF GRASHOF NUMBER 
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Figure 11-3: THERMAL CONDUCTIVITY OF VARIOUS GASES 
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Figure 11-4: CANDIDATE SCALE RATIOS - TEMPERATURE PRESERVATION 
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The model scale ratio and pressure ratio are fixed by the gas in the prototype 
and the gas chosen for the model* For example a prototype using air at 
atmospheric pressure could be simulated with a 3/8 scale model with Freon 11 
at a pressure of about 0.5 atm. 

Application of the "temperature preservation technique" to thermal scale model- 
ing of systems involving convective heat transfer must be evaluated separately 
for each particular system of interest. 

Practical problems that might arise in using this scaling technique include: 
gas availability, gas toxicity, compatibility between gas and model materials, 
liquefaction of the gas and lack of reliable thermophysical property data for 
the gas . 

The application of the temperature preservation scaling technique to the cabin 
atmosphere-cabin wall thermal interface is severely limited by the availability 
of suitable gases. Of the gases considered, only Freon 11, Krypton and Xenon 
allow scale ratios less than 0.5. Krypton and Xenon are rare gases and are not 
readily available. Freon 11 is available, however, it liquefies at 75°F at 
atmospheric pressure. Low pressure operation of a 3/8 scale model would be 
feasible using Freon 11 since its boiling point can be reduced to about -30°F 
by lowering the pressure to 1 psia. 

II. 4. 3 Scaling Compromises 

The problems associated with the "modified material preservation and "temperature 
preservation" scaling techniques limits their usefulness* Consequently other 
techniques must be sought for thermal scale modeling of systems involving 
convection. 

By allowing scaling compromises it may be possible to achieve adequate thermal 
similtude while preserving both gas and temperature in the scale model. Thermal 
similtude for the wall is achieved by preserving the convective heat transfer 
term in equation (5) . This term may be written as 
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Where q is the convective heat transfer rate per unit area of surface, 
conv 


The convective scaling criteria can then be written as 



(17) 


The necessary and sufficient conditons for satisfying equation (17) are, of 
course, given in Table II-l . However, two techniques which approximately 
satisfy equation (17) are suggested by rewriting this equation with the convec- 
tive heat transfer in terms of either heat transfer coefficient and the temper- 
ature difference between gas and wall, or mass flux and the gas temperature change 
through the system, i.e., either 



or 



(19) 


Temperature preservation using the same gas in model and prototype requires that 
both the heat transfer coefficient and the mass flux be preserved, i.e.. 




( 20 ) 


Even though the heat transfer coefficient and mass flux cannot be preserved 
simultaneously, preservation of either may result in adequate thermal similtude. 

Preservation of the heat transfer coefficient may be achieved if its dependence 
on the system parameters is known. In practical applications the heat transfer 
coefficient is usually based on analytical, semi-emperical or emperical results 
for simple flow fields (e.g. , flow over a flat plate). These results typically 
relate the Nusselt number to the Reynolds number for forced convection cases by 
an equation of the form 

Nu — C P& Pr (2i) 
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and to the Grashof number for free convection cases by an equation of the form 


A lu = c Gr n Pr M 


( 22 ) 


The values of the constants in these equations depend on the flow conditions 
(turbulent or laminar), boundary conditions and the system geometry. If the 
system to be modeled can be described by these relationships and if the constants 
are known, then the heat transfer coefficient can be preserved as follows: 


a. Laminar flow in model and prototype or turbulent flow in model and prototype 


Re#> * Re, 


Gr 


m 


(L* ) n 

\Lp) ' 


for forced convection and 


for free convection. 


b. Turbulent flow in prototype and laminar flow in model 

Re, 


m 


( Lp C u ) 


J2l ("Ll. zSHi-\ 

L Re n ^ Pr n -' ; 


forced convection and 


Gr m = §f)^Sr £ Pr ( 


nhr-jri L ) 

n L J 


for 


for 


free convection. 


c. Flow transition in either model or prototype 

If transitions between laminar and turbulent flow occur in either the model 
or prototype then it will not be possible to preserve the local heat transfer 
coefficient. 

Preservation of the heat transfer coefficient is most easily achieved for case 
(a) since only the exponents (n) are needed to determine the scaling criteria, 
whereas, case (b) also requires the constants (C) and exponents (m) to be known. 

The mass flux is easily preserved by adjusting the Reynolds number in the model 
such that 
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The free convection heat transfer is preserved for the mass flux as well as for 
the heat transfer coefficient preservation techniques* 


Preservation of either mass flux or heat transfer coefficient results in some 
loss of thermal similtude. If the mass flux is preserved then the ratio of 
forced convection heat transfer coefficient between model and prototype is 
given, for case (a) above, by 

** *= (l *; n " 

and conversely, if the forced convection heat transfer coefficient is preserved, 
then the ratio of mass fluxes between model and prototype is given by 

(f >&)*■-* ( l *)^ (25) 

Mass flux preservation tends to preserve the gas temperature change through the 
system, whereas, heat transfer coefficient preservation tends to preserve the 
temperature difference between the wall and gas. 

The degree of thermal similitude achieved with the mass flux and heat transfer 
coefficient preservation scaling techniques may be approximated using the 
following simplified analysis: 


Consider gas flow through an enclosure which is subject to uniform 
surface heating and coupled to the external environment through a 
uniform conductance per unit surface area. The thermal balance for 
the enclosure is given by 

QA S = KA $ (7~~ +JA s (T~Tf) (26) 

where 

Q = surface heating rate /unit area 

A = surface area 
s 

K = conductance per unit area to environment 
h = heat transfer coefficient 
T = enclosure temperature 
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T 

e 

T 

g 


= environment temperature 
- average gas temperature 


The average gas temperature may be related to the inlet gas temperature by 
equating the convective heat transfer rate to the energy transport rate 

/4(' r -'^) =vvc r^-r i 'j < 2? > 

where 


W = gas flow rate 
C = gas specific heat 


Setting, T = and solving for T gives 

g * # g 



t $*+ 


/ ■+ 


A s y- 

2 . »VC 




Using this result in equation (26) and solving for T yields 


" -7? = 


7* — Tr. 


/ +■/< 


j4s 

2ttC. 


.) 


where the reference temperature is the enclosure temperature 
when there is no gas flow and is given by 



(28) 


(29) 


(30) 


Assuming the free and forced convection heat transfer coefficients to 
be additive the heat transfer coefficient may be written as 



(31) 
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Using this definition of h and writing WC in terms of RePr equation 
(29) becomes 


T- T r — 


where 




/ + 


(%■)( 


C*T 




t 

Cfo”Pr m +C&r*'Pr 



(32) 


Ot = A /cross sectional flow area 
s 

The prototype and scale model temperatures can be written (for Pr = 1 and 
preservation of the free convection coefficient) as 


/ 


Prototype 



TW 

Scale Model (Mass Flux Preservation) 

TL- 


3 -f 


(KL P \r« 

(^/Ueep 


ce«J + c 




OC 


c ferv + c ' Gr e 




(33) 


(34) 


Scale Model (Heat Transfer Coefficient Preservation) 




ce«p +c,sr P 




-l 


(35) 


Figure II-5 shows the calculated temperature differences between model and 
prototype as a function of Reynolds number and Grashof number for the two 
scaling techniques . These calculated results used the following parameters 

m __ / 

L P ~ 5 

n " 0.5, n^ = 0.25 and 0^ = 0 = 0.33 (typical laminar flow values) 



(representative value for manned spacecraft 
cabin) 



(value for cylinder with length = radius) 
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The degree of thermal similitude achieved with either scaling technique is 
greatest at large Reynolds numbers where the wall temperature approaches the 
gas temperature and at small Reynolds numbers where the wall temperature 
approaches the reference temperature. The heat transfer preservation 
technique gives the best results for pure forced convection (Gr=0) and mass 
flux preservation gives the best results when free convection dominates. 


These scaling compromises also affect the transient response of the spacecraft 
cabin wall. The transient response may be calculated by substituting 
for T^ in equation (32) and solving for T, assuming T~T^ at t=0> 


T-T r 



(36) 


where 

^ “ (zj^) {jfePr * CRe”Pr'*+C ) <3r*'Pr m ' ) 

and H=thermal mass per unit area of wall 


The characteristic response time is then 



2 

Noting that (K/M) corresponds to the (k/pCL ) in the scaling criteria, 
equation (38) may be written in terms of nondimensional time as 



(39) 


The characteristic response time ratio between model and prototype for 
changes in convective heat transfer is then given by 


T 




(40) 
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This response time ratio for the 1/5 scale model spacecraft parameters 
is shown in Figure II-6 as functions of Reynolds and Grashof numbers for 
the two scaling techniques . 

As would be expected the relative model response is slow with heat transfer 
coefficient preservation and fast for mass flux preservation. 

The heat transfer coefficient preservation technique generally gives better 
transient simulation except for the lower Reynolds numbers when free convection 
effects dominate. 

The calculated results for this simple model indicate that the prototype temp- 
erature remains bounded by the scale model temperatures obtained using mass 
flux and heat transfer coefficient preservation. 

Whether mass flux preservation or heat transfer coefficient preservation gives 
better thermal similitude depends on the nature of the system to be modeled. 

Systems with mass flow rates large enough to give a small gas temperature change 
should be modeled well using the heat transfer coefficient preservation technique. 

If the free convection mode of heat transfer dominates then the mass flux 
preservation technique should give good results . 

II. 4. 4 Nusselt Number Preservation 

The Nusselt number preservation scaling technique would use the thermal scale 
model to experimentally determine the Nusselt number for the system as a 
function of Reynolds number and Grashof number. (If the model uses a different 
fluid the Prandtl number effects would also have to be determined) . This 
functional relationship would then be used in conjunction with a thermal math 
model to predict the prototype performance. 

The Nusselt number is preserved if thermal and dynamic similitude exist between 
the fluid elements of the systems. Thermal and dynamic similitude may be 
achieved by preserving Re and (Gr0^) . This keeps the nondimens ional velocity v^ 
and fluid temperature invariant for similar systems. (see equations 10 and 

11). The Reynolds number is easily preserved in the scale model, however 
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preservation of the term depends on the magnitude of the convective 

heat transfer. In the limit of no convective heat transfer — | 

and in the limit of only convective heat transfer && * L 4 * 

Consequently the model pressure required to preserve the GtQ^ term (i.e, the 
free convection effects) is within the limits 

(l!*)"* <«i) 

The use of Nusse.lt number preservation scaling technique can be summarized as 
follows: 

Use temperature preservation technique for radiation - conduction 
aspects of thermal scale model. 

° Develop math model for radiation-conduction aspects of scale model 

° Test model without convection and use results to upgrade math model 

° Test model with convection at various pressures for each Reynolds 

number . 

° Use the experimental results in conjunction with math model to deter- 
mine the convective heat transfer. 

° Develop Nusselt number correlation with Reynolds and Grashof numbers 

° Expand math model to include convection effects and use in conjunction 
with the Nusselt number correlation to predict prototype performance. 

The Nusselt number preservation technique avoids the problems inherent in the 
other scaling techniques, i.e., high temperatures and heat fluxes required for 
the "modified material preservation" , the suitable gas required for temperature 
preservation", and the uncertainties of "scaling compromises." 

However thermal scale modeling using Nusselt number preservation requires the use 
of experimental research techniques whereas use of the other scaling techniques 
requires only fabrication and testing of the scale model. 
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11*4.5 Simulation of Zero Gravity Conditions 


In order to simulate the pure forced convection which occurs in a zero gravity 
field, free convection must be reduced to a negligible amount in the ground test 
simulation. Free convection effects on forced flow are generally related to the 


Gt 

magnitude of for flow over a vertical surface. 


vection dominates and if 


Gf 

Be 4 


If — 2 ^ 1 forced con- 


»1 free convection dominates. 


The Grashof number to Reynolds number squared ratio may be written as 

— = gejrCszJhl ( 42 , 

Using a characteristic length of 10 feet, temperature differences between wall 
and gas of 1. to 10. °F and gas velocities of 15 to 45 ft/min as representative 
for manned spacecraft yields values for ° ^ to WO for g roun d tests 

of manned spacecraft configurations. This indicates that free convection effects 
are important, and in fact may dominate, in ground tests of prototype manned 
spacecraft. 


The reduction of free convection effects in thermal scale model tests depends on 
the scaling technique, scale ratio and gas density. At a given gas pressure the 
value of — * decreases with smaller scale ratios for the Modified Material 
Preservation, Nusselt Number Preservation and Mass Flux Preservation Scaling 
techniques and increases for the Heat Transfer Coefficient Preservation Technique. 


Gp 

The value of can t> e reduced by decreasing the gas density while preserving 

the Reynolds number. This means that the gas velocity must be increased to 
compensate for the reduced density. The maximum allowable gas velocity is 
limited by the onset of compressible flow effects. This reduction in “ffet 
could also be obtained in the prototype tests. 


Car 

The scale model to prototype — ^ ratio may be written as 


(&/& )*= '* ( L *) '(***) W '*(' T *)'\ p *) 


(43) 
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This may be written for the various scaling techniques as follows: 


// 


Modified Material Preservation 


Temperature Preservation 




Mass Flux Preservation 


L *(p*y 


Heat Transfer Coefficient Preservation (Laminar Flow) 


C Gr /Re z T = ( L K )~'(P“) 


(44) 


(45) 


(46) 


(47) 


Nusselt Number Preservation 




(48) 


The free convection effects in manned spacecraft scale model tests should be 


eliminated by reducing 


Sr 


by a factor of 1000. This corresponds to 


values of 0.001 to 0.1 for the scale model tests. The pressure ratio required 
for this reduction in free convection effects is shown in Figure II-7 for the 
various scaling techniques . Even though these required pressure levels result 
in cabin gas velocities well within the incompressible flow regime, problems in 
sizing the gas supply and return lines may arise for the techniques requiring the 
lower pressure levels. 


II. 4. 6 Manned Spacecraft Applications 


The "modified material preservation 11 scaling technique applied to manned space- 
craft is limited, by the increased model temperature, to scale ratios greater 
than about 0.5. The increased model pressure requirements preclude simulation of 
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1-g prototype gravity field. This technique is well suited for zero gravity 
simulation. 

The "temperature preservation" scaling technique is limited by availability of 
"suitable gases." A 3/8 thermal scale model of a manned spacecraft appears to 
be possible using Freon 11. Simulation of a 1-g prototype gravity field is 
marginal because of possible liquification of the Freon gas. The low model 
pressure required to simulate zero gravity conditions may result in design 
problems for the gas supply and return lines. 

Scaling compromises using the mass flux and/or heat transfer coefficient 
preservation scaling techniques may give adequate thermal similitude without 
limiting the model scale ratio. Use of both techniques should bound the pro- 
totype spacecraft performance. Both techniques are suited for simulation of 
1-g prototype conditions. However the low model pressure required for simula- 
tion of zero gravity conditions using "heat transfer coefficient preservation" 
may result in design problems for the gas supply and return lines. 

The "Nusselt number preservation" scaling technique is unsuited for simulating 
1-g conditions in the prototype spacecraft due to the increased pressure 
requirement. This technique is well suited for simulating zero gravity condi- 
tions without limiting the model scale ratio. 

Since "scaling compromises," using both mass flux and heat transfer coefficient 
preservation, and Nusselt number preservation allow a greater range of thermal 
scale modeling applications these techniques were chosen for an experimental 
investigation of thermal scale modeling applied to a radiation-conduction- 
convection system. 
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II. 5 EXPERIMENTAL INVESTIGATION 

11. 5.1 Objective 

The objective of the investigation was to experimentally demonstrate the applica- 
tion of the selected scaling techniques to thermal scale modeling of systems 
involving radiation-conduction and convection, in particular, to systems represen 
tative of the manned spacecraft "Cabin-Atmosphere - Cabin Wall Thermal Interface. 

In order to accomplish this objective the following plan was developed: 

° Select a model configuration in which radiation, conduction and 
convection are important* 

° Construct full scale and 1/4 scale models using temperature preser- 
vation scaling criteria for the radiation and conduction aspects of the 
configuration. 

° Develop a Thermal Math Model (TMM) for the configuration. 

° Test the 1/4 scale model without convection at various heating rates 
and upgrade and/or verify the radiation-conduction aspects of the TMM. 

° Test the models for various heating rates and gas convection conditions 

° Use verified TMM in conjunction with 1/4 scale model test results to 
evaluate the Nusselt number preservation scaling technique. 

° Correlate the 1/4 scale model and full scale model test results to 
evaluate the mass flux and heat transfer coefficient preservation 
scaling techniques* 

11. 5. 2 Model Configuration 

The model configuration chosen for the experimental investigation is shown 
in Figure II-8. The model basically consists of two concentric cylinders with 
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end plates closing off the annular region between the cylinders. Thermal grad- 
ients are generated by heating the center section of the inner cylinder and 
cooling the center of the outer cylinder. 

The heating is provided by an electrical heater mounted on the inside of the 
inner cylinder and the cooling by a cooled fin attached to the center of the 
outer cylinder. 

The exposed surfaces in the annular region are painted black to enhance the 
radiation heat transfer and to provide the same radiation properties for full 
and 1/4 scale models. 

Forced convection is provided by air flow in the annular region. Air is supplied 
to and exhausted from the model through air distribution plugs in the air inlet 
and outlet sections of the inner cylinder. 


This model configuration was chosen for the following reasons: 

° Radiation, conduction and convection modes of heat transfer are 
important in the thermal response. 

° Symmetry of model minimizes the amount of instrumentation 

and the size of the TMM and also allows the asymmetric convective 
heat transfer effects to become apparent. 

° Air flow distribution is similar to that for manned spacecraft 
configurations . 

° Expensive space chamber testing is avoided by insulating the model 
and testing at ambient conditions. 

The annular region can be evacuated for radiation-conduction tests 
and can be pressurized to increase the free convection effects. 

° The annular flow region can be kept free of instrumentation and heater 
leads . 
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o The inner cylinder can be evacuated to reduce the extraneous heat 
transfer in this region. 

The basic structure of the full scale model uses 6061-T6 aluminum and measures 
approximately 2 feet in length and diameter. The 1/4 scale model uses type 304 
stainless steel and measures approximately 6x6 inches. The basic length 
dimensions for the two models are shown in Figure II-9. The inner cylinders 
have thicker walls at the center and end sections , where the heater and air 
distribution plugs are installed. The model material thicknesses and the scaling 
parameter for two dimensional conduction are given in Table II-3. The conduction 
scaling parameter agreement is generally within 10 percent over the temperature 
range of 30 to 350°F. These model dimensions reflect the requirements of 
pressurization and evacuation while maintaining relatively thin walls to keep 
conduction from dominating the heat transfer processes. 

II. 5 .3 Basic TMM and Preliminary Analysis 

A TMM was developed for the basic configuration and preliminary analyses were 
made to verify the thermal similtude between the full scale and 1/4 scale models 
and to determine the experimental accuracy required to implement the Nusselt 
number preservation scaling technique. 

II. 5. 3.1 Basic TMM 

The basic TMM nodal network is shown in Figure 11-10. There are 12 equally 
spaced nodes (401-412) on the outer cylinder and 12 (301-312) on the inner 
cylinder. Both end plates have 4 equally spaced nodes (101-104 on the air out- 
let end and 201-204 on the air inlet end) . The fin node is 413 and the boundary 
node at the outer edge of the fin is node 414. The inlet gas temperature 
boundary node is numbered 500 and the outlet gas node is number 513. The 
remaining gas nodes (501-512) are equally spaced axially in the annular region 
between the inner and outer cylinders . 

The radiation exchange factors (script F’s) between the interior solid nodes 
(101-104, 201-204, 301-312 and 401-412) were calculated using Boeing’s General- 
ized Radiative Heat Transfer Program (AS2814) (Reference 8) for the surfaces 
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TABLE II- 3 MODEL MATERIAL THICKNESSES AND SCALING PARAMETERS 


Component 

Model 

Thickness 


Scaling Parameter* 




d, inches 


BTU/HR- 

°F 





30 °F 

75°F 

200 °F 

350°F 

Outer Cylinder 

Full Scale 

0.063 

.502 

.508 

.521 

.527 


1/4 Scale 

0.040 

. 425 

.446 

.496 

.531 

Inner Cylinder 

Full Scale 

0.050 

, 399 

.404 

.414 

.418 

Thin Section 

1/4 Scale 

0.032 

.339 

.366 

.396 

.424 

Thick Section 

Full Scale 

0.125 

.998 

1.010 

1.035 

1.046 


1/4 Scale 

0.085 

.900 

.946 

1.052 

1.128 

End Plates 

Full Scale 

0.375 

3.00 

3.02 

3.10 

3.14 


1/4 Scale 

0,250 

2.65 

2.78 

3.10 

3.32 

Cooling Fin 

Full Scale 

0.071 

2,51 

2.52 

2.54 

— 


1/4 Scale 

0.063 

2.44 

2.50 

2.56 


^Scaling parameter 

2 

is kd/N for wall 

components , where 

k ~ thermal conductivity, 


d = wall thickness 

, N = scale ratio. 

For the cooling 

fin the 

scaling parameter 

is 




where R^, = outer radius of fin, (R + d Q ) = outer radius of outer cylinder. 
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HEAT SOURCE AT NODES 306 & 307 
NODES 414 & 500 ARE BOUNDARY NODES 
O SOLID NODES 
• GAS NODES 


Figure U-10: BASIC THERMAL MA TH MODEL 
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painted with black paint 0.88) . These calculations resulted in 496 
SAF (aAF) values which were input into the BETA (Boeing Engineering Thermal 
Analyzer) program (Reference 9) * 

The gas flow was treated as a source at the gas nodes. The value of the source 
at node i is given by = wc (T^^ - T^) where w * mass flow rate and c - 

specific heat of the gas, e.g., = wc (T*.^ “ T^^) , conductors 

were input as a function of temperature and heat sources input at nodes 306 and 
307. 

The conductors (hA) between the wall and the gas were based on an approximation 
of the combined free and forced convection and the gas conduction heat transfer. 
The forced convection heat transfer coefficient was assumed to be that given for 
laminar flow over a flat plate 

4 (x) = 0.332(A) *4 Pr*(?Y 

and free convection heat transfer assumed to be 

Nu - 0 . 04 £r^ 

or basing the Grashof number on the annular spacing (R q - R^) and the local gas 
to wall temperature difference 



The total heat transfer coefficient used in the analysis consists of the sum 
of the forced convection, free convection and gas conduction 
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The + sign is used when the buoyancy acts with the flow and the - sign when 
it acts against the flow* Along the end plates the buoyancy term is not 
included and the forced convection term is modified to account for the radial 
flow of the gas along these plates. These convective heat transfer relationships 
are only a reasonable guess as to how the convective heat transfer occurs, 
however; they serve their purpose in the assessment of the experimental 
configuration. 

II 5.3.2 Preliminary Analyses 

Thermal analyses, using the math model were made for the full scale model and 
the 1/4 scale model for the following conditions: 


Full 

Size Model (heat input 

1000 Btu/hr) 

1 . 

Evaculated (radiation ■ 

- solid conduction) 

2. 

Low Pressure (#1 + gas 

conduction) 

3. 

1 atm Pressure (#2 + free convection) 

4. 

Gas Flow (#3 + forced < 

convection) 


a. Re = 1180 

V = 6 . f pm 


b. Re = 4720 

V = 24. fpm 


c. Re = 18,900 

V = 96. fpm 

1/4 

Scale Model (heat input 

= 62.5 Btu/hr) 

1 . 

Evacuated 


2. 

Low Pressure 


3. 

1 atm Pressure 


4. 

Gas Flow 



a. Re = 1180 

V - 48. fpm 


b. Re - 4720 

V - 192. fpm 


c. Re = 18,900 

V = 788. fpm 


5. #3 and #4 repeated at pressures of 4. , 8. and 16. atmospheres. 

These cases were analyzed for air inlet (where applicable) temperatures of 
35 and 80°F. All cases used a cooling fin temperature of 35°F. 
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The calculated temperature distributions for an air inlet temperature of 80 °F 
are shown in Figures 11-11 and 11-12 for the full scale and 1/4 scale models 
respectively . The full scale model temperatures are somewhat lower than those 
of the 1/4 scale model for the radiation-solid conduction case. It was originally 
planned to use 6061-T4 alloy aluminum for the full scale model, however, it turned 
out that 6061-T6 alloy was used instead. This switch in full scale model alloy 
resulted in model conductances about 10 percent too low. Figure 11-13 shows the 
calculated temperature differences between the 1/4 scale and full scale model for 
both 6061-T4 and 6061-T6 aluminum alloys. The good agreement with the 6061-T4 curve 
reflects the intended use of this alloy. The calculated 1/4 scale model temperatures 
average about 7.5°F warmer than the full scale model constructed of 6061-T6 
alloy. The 1/4 scale model fin design was subsequently modified to properly scale 
the full size model. This change lowers the model temperature by about 5°F. 
Consequently the imperfect thermal conduction scaling is expected to result in a 
1/4 scale model temperature that is too high by about 2.5°F. This represents the 
largest differences to be expected since the cases analyzed are for the highest 
heating rates anticipated during radiation-conduction testing. 

A subroutine to the TMM BETA program was written to evaluate the experimental 
accuracy required to implement the Nusselt number preservation scaling technique. 

This subroutine calculates the convective heat transfer coefficient at each node 
using a given temperature distribution. When the calculated temperature distri- 
bution is used the input heat transfer coefficients are calculated. By using the 
calculated temperature distribution with the temperatures rounded off to varying 
degrees of accuracy, the effect of temperature measurement accuracy on the 
calculated heat transfer coefficients is determined. 

The effect of temperature measurement accuracy on the heat transfer coefficient 
determination was investigated by rounding off the calculated temperatures to 
the nearest 0.1, 0.2 and 0.4 degrees and using the subroutine to calculate the 
heat transfer coefficient. Figure 11-14 shows a typical comparison between the 
input heat transfer coefficients and those calculated for the various tempera- 
ture accuracies. The accuracy of the temperature measurements is most critical 
in regions where the thermal gradient is small and/or where the wall temperature 
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is near to the gas temperature. The results shown in Figure 11-14 indicate 
that the heat transfer coefficient may be determined to within 10 to 15 per- 
cent with temperatures measured to +0 .1°F, to within about 30% for +0 . 2° F and 
to within about 80% for +0.4°F. These results are for an inlet gas temperature 
of 35° F. The results for an inlet gas temperature of 80° F show larger uncer- 
tainties in the determination of the heat transfer coefficient for the outer 
cylinder regions where the gas temperature and wall temperature are nearly equal. 
Consequently it was decided that a temperature measurement accuracy of at least 
+0.1° F is required to implement the Nusselt number preservation scaling 
technique. 

11.5.4 Thermocouple Calibration 

Forty-nine chromel-cons tantan thermocouples were calibrated for the full scale 
model and a like number for the one-quarter scale model. The wires are AWG #36 
(0.005 in. diameter) with a double layer of stranded fiberglas insulation. 

All the wires were taken from the same spool. The junctions were made by 
twisting the leads together and silver brazing. 

The calibration consisted of an absolute calibration and a relative calibration. 
The absolute calibrations were conducted, on samples from the spool, at the 
Boeing Metrology Laboratory. This calibration has NBS traceability and claims 
an accuracy of 4p*04°F for any particular calibration point. This calibration 
data is presented in Table 11-4... The standard EMF is based on tabulated results 
of reference 10. The relative calibration consisted of fixing the forty-nine 
thermocouples earmarked for a particular model on a copper slug. The slug was 
then put in an insulated oven fixture and all the wires were checked against 
each other over the entire test temperature range. Other than the transferring 
of the test junctions from the copper slug to the model, all the thermocouple 
circuitry for the relative calibration was kept for the model tests. The oven 
and the thermocouple circuitry are shown in Figure 11-15. 

Since the preliminary analyses showed that a temperature roundoff (simulated 
thermocouple error) larger than 0.1°F began to cause scatter in the calculated 
convective heat transfer coefficient, it was decided to correct all thermocouples 
to at least this accuracy. For chromel-cons tantan wire at the test temperatures. 
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TABLE II-4 THERMOCOUPLE ABSOLUTE CALIBRATION 


Thermocouple EMF (above Ice point) 


Temperature 

Standard 

Measured 

Correction 

(°F) 

(m V) 

(m V) 

(Ji V) 

32 

0 

0 

0 

100 

2.2753 

2.2708 

+4.5 

200 

5.8724 

5.8598 

+12.6 

300 

9 . 7112 

9 .6950 

+16.2 

400 

13.7518 

13.7378 

+14.0 
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Figure 11-15: THERMOCOUPLE CALIBRATION SET-UP AND CIRCUITRY 
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0 * 1°F is approximately equal to 3.5 pV of thermocouple emf. This allowed 
grouping the thermocouples for purposes of generating correction curves. For 
the full scale model the grouping was six groups of 6, 6, 12, 13, 6 
and 6. For the quarter scale the grouping was seven groups of seven each. Each 
group can then be represented by its own mean and corrected to the mean of the 
forty-nine wires. The absolute calibration can then be added to the mean of the 
forty-nine wires. In this way each wire is accurately related to every other 
wire in the model (+3 yV) and is also accurate to within approximately 0.5°F 
of absolute temperature. Note that the +3 yV is the maximum error between any 
two thermocouples and represents a simulated round off of individual temperatures 
in the math model of approximately 0.05°F. The relative calibration curves for 
the full scale and the one-quarter scale model are given in Figure 11-16 . 

II. 5. 5 Model Design, Fabrication and Instrumentation 

The basic model designs (outer cylinder, cooling fin, end plates and inner 
cylinder) were outlined in Section II. 5. 2. The outer cylinders were fabricated 
from sheet material which was rolled and seam welded to form the cylinders. The 
cooling fins were also made from sheet material. The horizontal fins were 
cut from a sheet and welded to the vertical fins which had been formed by rolling 
and seam welding. These fin structures were then welded onto the centers of the 
outer cylinders. The end plates were cut from plate material. The inner 
cylinder for the 1/4 scale model was machined from a 2.0 inch O.D. by 0.109 inch 
wall tube. The inner cylinder for the full scale model was fabricated using 
tubing and sheet material. The thicker sections of the inner cylinder were cut 
lengths of 8.0 inch O.D. by 0.125 inch wall tubing and the thinner sections 
were formed from 0.050 sheet material by rolling and seam welding. These sections 
were then welded together to form the inner cylinder. 

Twelve equally spaced holes were cut at each end of the inner cylinders to 
provide air inlet and outlet passages. These holes were cut such that one side 
would be flush with the end plates in the final model assemblies. The holes are 
slots with semicircular ends. The full scale model holes measure 1.344 inches 
long by 0.625 inches high and the 1/4 scale model holes 0.336 inches long by 
0.156 inches high. 
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The air distribution plugs were designed to provide air to the annular region 
of the models. These plugs also seal off the inner cylinder while providing 
instrumentation and heater lead feed throughs, as well as providing the capa- 
bility for evaluating the inner cylinder. No attempt was made to thermally 
scale model these air distribution plugs since their effects on the thermal 
performance of the models should be minor. The flow passage geometry was scaled 
to preserve dynamic similarity. Both full scale and 1/4 scale model plugs 
were made of phenolic material. Figure 11-17 shows the 1/4 scale model air 
distribution plug that has the heater and thermocouple feed through passages. 

The other plug is identical with the exception of having only two pass through 
holes 180° apart which are counterbored to accept 1/4 inch O.D. vacuum lines. 
Twelve equally spaced 0.156 inch diameter holes drilled radially provide the 
air distribution passages. The air supply (or exhaust) is provided by a 
hole in the center of the plug which is threaded to accept a standard fitting 
for 1/2 inch O.D. tubing. The full scale model air distribution plugs are 
similar to those for the 1/4 scale model except that only two feed through? 
were required for the instrumentation and heater leads. 

The model heaters were designed to provide a uniform heat flux to the center 
section of the inner cylinders. The full scale model heater was fabricated 
by uniformly wrapping 0.001 x 0.062 inch Nichrome V ribbon on a flat 0.015 inch 
thick phenolic sheet which was then bonded to a 0.125 inch silicone rubber 
backing sheet. The exposed heater ribbon was then coated with a layer of RTV 
silicone rubber approximately 0.006 inches thick. The heater assembly was 
sized such that when rolled into a cylinder it would fit the inside of the center 
section of the inner cylinder. The 1/4 scale model heater was fabricated by 
wrapping 28 turns of Nichrome wire around a threaded heater core that had been 
machined from a "Transite” block. A 0.059 inch thick layer of RTV silicone 
rubber was then bonded on the outer surface of the assembly. This assembly was 
designed to be a press fit inside the center section of the inner cylinder. 

Prior to installing the heaters, 24 holes (0.015 inch diameter) were drilled in 
each inner cylinder for the thermocouple instrumentation. Two holes were 
drilled 180° apart at 12 axial locations along the inner cylinders. These 
axial locations correspond to the center of each inner cylinder node of the TMM. 
An exception to this positioning occurs for the air inlet and outlet sections 
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where the air holes interfere* The thermocouple holes for these sections are 
75 percent of the node length away from the end plates instead of being at the 
node center. 

The four thermocouples for the center section of the inner cylinders were 
installed before the heaters were placed in position. These thermocouples (as 
were all of those for the inner cylinders) were installed by inserting the 
twisted and silver brazed junctions in the holes, securing them in position by 
swaging, filing off the junction excess protruding beyond the outer surface of 
the cylinder and sealing the holes with epoxy resin. 

Figure 11-18 shows a photograph of the 1/4 scale model components prior to 
installation of the heater inside the inner cylinder. This figure shows the 
two end plates, the outer cylinder assembly with the copper coiling coils brazed 
to the fin, the inner cylinder with the four center section thermocouples 
installed, the air distribution plugs and the heater assembly. The heater 
assembly is shown without the 0.059 inch layer of RTV silicone rubber. This 
layer was bonded to the heater in two pieces leaving two gaps to clear the center 
section thermocouple leads* After the heater was pressed into place, these gaps 
were filled with RTV using a hypodermic syringe. The full scale model heater was 
installed by coating the heater surface with RTV, inserting the heater in 
place inside the center section of the inner cylinder and holding it tightly in 
position against the cylinder wall, using an inflatable bladder, until the RTV 
cured. After the RTV cured the gap in the heater assembly was filled in with 
RTV. The full scale model heater design and installation is shown in Figure 
11-19. Following the heater installation, the remaining inner cylinder thermo- 
couples were installed and "bench tests" were made to check out the thermocouple 
and heater installations. 

One end plate was welded to the inner and outer cylinder assemblies, and the 
radiation surfaces were cleaned and painted with two coats of Fullers Metal 
Etch Primer (3811 black with 3816 primer). Figure 11-20 shows the full scale 
model at this stage of assembly* The emittance data for the black paint was 
measured for various painted 6061-T6 aluminum and 304 stainless steel samples. 
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Figure 11-19: FULL SCALE MODEL HEATER INSTALLATION 
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FIGURE 11-20 FULL SCALE MODEL BEFORE FINAL WELDING 
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These data are presented in Table II-5. The preliminary analyses assumed an 
emissivity of 0.88. The welding of the remaining end plates to the inner and 
outer cylinders completed the basic structure assemblies. The "0" ring seal 
surfaces at each end of the inner cylinders were machined after the model 
assemblies were welded. 

The inner cylinder assemblies were completed with the installation of the air 
distribution plugs. Figure 11-21 shows the inner cylinder assembly of the 1/4 
scale model. The air distribution plugs are held in position by a 0.75 inch 
diameter phenolic rod which keeps external pressure from pushing the plugs 
further into the cylinder. The full scale model uses a 2 inch O.D. by 1/16 inch 
wall phenolic tube to support the air distribution plugs. At this point, a 
recheck of the thermocouples revealed that many of the 1/4 scale model thermo- 
couples were broken. Because of the limited working space inside the inner 
cylinder, the 1/4 scale model was disassembled and reins trumen ted . The welds 
between the end plates and outer cylinder were cut and the inner cylinder assembly 
with attached end plates was removed from the outer cylinder assembly. The heater 
and thermocouples were removed, the old thermocouple holes welded closed, and new 
holes drilled. The thermocouple leads were shortened to reduce the possibility 
of subsequent damage during handling of the model and new thermocouple junctions 
fabricated and installed as shown in Figure 11-22. Since the old thermocouple 
leads were used for the new junctions the effect on the calibration curves is 
insignificant. After brazing the thermocouple junctions the bared lead wires were 
encapsulated with RTV (silicone rubber) to form tabs about 3/8 inch square by 
0.05 inch thick with the brazed functions emerging normal to the tabs. The 
junctions and adjacent tab surfaces were coated with RTV and the junction were 
inserted into the drilled holes of the inner cylinder and held in place. The 
curing of the RTV then bonded the thermocouple assembly in place. The thermocouple 
leads were wrapped a distance around the inner cylinder wall in the plane of the 
junction (to reduce the heat leak from the junction) and then routed through the 
pass through holes in one air distribution plug. The thermocouple and heater leads 
were sealed in the pass through holes with "Scotchcast 5. 11 Prior to the 
installation of the air distribution plugs, the inner cylinder was filled with 
"Micro Quartz" insulation to reduce heat transfer inside the inner cylinder. 

Before rewelding the end plates to the outer cylinder an extensive vaccum leak 
check was made for the inner cylinder using a helium leak detector. However, 

even after sealing all of the leaks that could be found, the minimum pressure 

-3 

that could be obtained in the inner cylinder was about 5 x 10 torr at room 

_3 

temperature and 80 x 10 torr with the heater at 300°F. 
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TABLE I 1-5 EMITTANCE DATA FOR PAINTED SAMPLES 
FULLERS METAL ETCH PRIMER (3811 BLACK WITH 3816 PRIMER) 


Sample No. 

Substrate Condition 

Number of 
Coats ( Sprayed) 

e 

e * ** 

1 

6061-T6, base 

1 

0.78 

0.73 

2 

6061-T6, base 

2 

0.88 

0.83 

3 

606 1-T6 , roughened 

1 

0 . 79 

0.75 

4 

6061-T6 , roughened 

2 

0.90 

0.85 

9 

304 CRES, base 

1 

0.78 

0.78 

10 

304 CRES, base 

2 

0.85 

0.88 

11 

304 CRES, roughened 

1 

0.75 

0.76 

12 

304 CRES, roughened 

2 

0.84 

0.88 


*Measured with a Lions Emissometer. 

** Value after baking cycle of: 

3 1/2 hours at 220°F 
I 1/4 hours at 350°F 

4 hours at 400°F 
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Figure 11-21: 1/4 SCALE MODEL INNER C YLINDER ASSEMBL Y 
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Figure 11-22: THERMOCOUPLE INSTALLA Tl ON FOR 1/4 SCALE MODEL INNER CYLINDER 
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The outgasing of the transit heater core and the RTV along with the small 
pumping line apparently limits the evacuation of the inner cylinder. 

Twenty-three thermocouples were attached to the outside of each model at 
positions corresponding to the TMM node locations. The junctions were spot 
welded to the model except those on the 1/4 scale model end plates were staked 
into drilled holes and those on the copper coiling tubes were brazed. The 
thermocouple installation corresponding to nodes 201-204 and 401-406 is shown 
in Figure 11-23 for the full scale model and Figure 11-24 for the 1/4 scale 
model. The thermocouple installation corresponding to nodes 101-104 and 407-412 
is shown in Figure 11-25 for the full scale model and Figure 11-26 for the 1/4 
scale model. The thermocouple installation for node 413 can be seen in Figures 
11-25 and 11-24 for full scale and 1/4 scale models respectively. Two thermo- 
couples were used for node 414, one at the inlet and the other at the outlet 
of the cooling coil. These may be seen in Figure 11-27 for the full scale 
model and Figure 11-26 for the 1/4 scale model. The thermocouple leads were 
routed around the models to minimize the heat leak from the junctions and were 
grouped into a single bundle of leads routed away from the model. The routing 
of the lead bundle away from the models can be seen in Figures 11-27 and 11-28. 
The correspondence between the thermocouple numbers and the TMM node numbers is 
shown in Figure 11-29 for the full scale model and Figure 11-30 for the 1/4 
scale model. 

Prior to testing, the models were insulated with polyurethane foam. The models 
were centered in cylindrical molds and the insulation was foamed between the 
models and the molds. The foam insulation is shown in Figures 11-31 and 11-32 
for the full scale and 1/4 scale models. The foam thickness was sized to 
preserve the one-dimensional heat transfer through the Insulation. Figure 11-33 
shows a photograph of the insulated full scale model. 

II. 5. 6 Model Tests 

The full scale model was tested at atmospheric pressure under free convection 
and forced convection conditions. Five heating rates were used in the free 
convection (no flow) tests. The forced convection tests used three heating 
rates and three flow rates. The flow rate tests were repeated for the highest 
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FIGURE 11-24 INSTRUMENTED Va SCALE MODEL 


360 









: : ■ 


A? 

' > i ,Sv« ' A i f \ " Av 

: ;, 


\ ; : 




, . 

A' ■■.■ ■ : 


: 


' ■ : 

mstm 


AA-A 


; ; : 

- 'i •• 


Mm 


mm 


'MmM'i'l 









o:!-; 




siliill 


wmmi 

r • V-vv 

S;«> 


. .: &§J 
j 


FIGURE 11-27 INSTRUMENTED FULL SCALE MODEL 
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1 /4-SECTION 



Figure 11-32: 1/4 SCALE MODEL INSULA TION 
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heating rate case with the flow direction through the model reversed to see 
if the flow direction has any effect on the data correlation. To minimize 
the number of 1/4 scale model tests the full scale model flow rates were 
chosen to differ by factors of four, i.e. , letting w = highest flow rate, 
w /4 = intermediate flow rate and w /16 = lowest flow rate. 

The 1/4 scale model tests consisted of radiation-conduction, free convection and 
forced convection tests. The purpose of radiation-conduction tests was to obtain 
data for upgrading the thermal math model. The free convection tests were at 
five heating rates corresponding to the full scale model heating rates with the 
heating rate scaled such that 

Qm — 'TtZ Qp 

These 1/4 scale model tests were made at pressures of 1/2, 1,2,4 and 8 atmos- 
heres . The heat transfer coefficient preservation scaling technique uses the 
1/2 atmosphere tests to simulate the full scale model tests. This assumes that 
the free convection Nusselt number is proprotional to the fourth root of the 
Grashof number. Based on this relationship the 1/4 scale model pressure required 
to preserve the free convection heat transfer coefficient is 



The Nusselt number preservation scaling technique uses all of the tests to 
determine the Nusselt number as a function of Grashof number. 

The 1/4 scale model forced convection tests were made to obtain data to be used 
for the Nusselt number, mass flux and heat transfer coefficient preservation 
scaling techniques. The forced convection Nusselt number was assumed to be 
given by flat plate laminar flow theory (i.e., Nusselt number proportional to 
the square root of the Reynolds number). Based on the relationship the 1/4 
scale model mass flow rate required to preserve the heat transfer coefficient 
is 


64- o 
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The mass flow rate required for mass flux preservation is 

= TZ % 

and that required for the Nusselt number preservation technique is 

Consequently three 1/4 scale model flow rates are required for each full scale 
model flow rate, however, since the full scale model flow rates differ by 
factors of four, only five different flow rates were required for the 1/4 scale 
model tests (i.e., w / 4, w /14, w /64, w /128 and w /512). The heat transfer 
coefficient and mass flux preservation tests were made at 1/2 atmosphere 
pressure and at the three heating rates corresponding to the full scale tests. 
The Nusselt number preservation tests were generally made at pressures of 

I, 2, 4, and 8 atmosphere pressure. The heating rates for these tests were 
not restricted to simulating the full scale model heading rates since this is 
not required for determining the Nusselt number as a function of Grashof number 
at given Reynolds numbers. 

A matrix of the model test runs is given in Table II-6. These are the test 
runs basic to the data correlation. The runs not shown are the various pre- 
liminary and checkout tests. 

II. 5. 6.1 1/4 Scale Model Tests 

II. 5. 6. 1.1 Preliminary Tests 

Prior to insulating the 1/4 scale model with polyurethane foam a checkout was 
made of the test setup and model instrumentation. A photograph of the test 
sqtup is shown in Figure 11-34. Figure 11-35 gives a schematic of the test 
setup. The 1/4 scale model is mounted on top of a vacuum cart which provides 
the vacuum pumping facility for the model. The heater power is provided by 
a D.C. power supply. The thermocouple readout apparatus is the same as that 
used for the thermocouple calibration. The cooling for the fin is provided 
by water which is colled in an ice-water bath. The ice-water bath is main- 
tained at the ice point by means of an ethylene glycol-water loop connected 
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**Heating rate given for 1/4 scale model. Heating rate for full scale model is 16 times the values given 
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Figure 11-35: PRELIMINARY TEST SETUP SCHEMA TIC 
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to a Freon refrigerator- Figure 11-36 shows another view of the test setup 
with the full size model along side the 1/4 scale model to give a direct 
size comparison. 

In order to checkout the test setup and model instrumentation, the 1/4 scale 
model was temporarily insulated with approximately one inch thickness of foam 
rubber, and preliminary tests were made with the vacuum system, cooling system 
and the heater turned on. These preliminary tests resulted in the following 
observations : 

o The out gassing of the paint in the annular region between the inner and 

outer cylinders slows the pump out of this region. 

o The cooling system is adequate to maintain the fin temperature at a 

constant value between 33 and 34°F. 

o Thermocouple numbers 14, 18 and 24 not giving any reading. Thermocouple 
number 12 reading erroneously. All other thermocouples give repeatable 
readings with the differences between adjacent thermocouples repeatable to 
within about 1 microvolt («*0.03°F). The heater section thermocouples 
show some discrepancies. Thermocouple number 6 is apparently reading 
correctly, however, numbers 7 and 19 appear to be reading a temperature 
away from the cylinder wall and closer to the heater core. 

o The copper vacuum lines on one of the air distribution plugs provide a 
significant heat leak path and need to be shortened. 

After these preliminary tests the air distribution plug containing the copper 
vacuum lines was removed and the tubes were shortened. With the air distribution 
plug removed an inspection of thermocouple numbers 12 and 24 revealed nothing 
abnormal with 12, however, 24 was broken near the junction (the junction x^as 
repaired and the thermocouple reinstalled) . 
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After insulating the model with polyurethane foam (as shown in Figure 11-32) 
it was re-installed in the test set up. The vacuum system was turned on and 
the model evacuated to about 0.2 p of Hg . in the region between the inner and 
outer cylinders. However, when the heater was turned on the pressure increased 
to about 1 p of Hg. due to the outgassing of the painted surfaces. In order to 
reduce the outgassing the heater was left on with the cooling water off in an 
attempt to ,f bake out" the system. After several days of "bake out" the pressure 
appeared to stabilize at 0.85 p of Hg. with the heater section at about 200°F 
and the outer wall at about 150°F. The heater was then turned off and the cooling 
water turned on to find the minimum attainable pressure with the model at below 
room temperature. The model had cooled to about 45°F and the pressure had dropped 
to 0.3 p of Hg. when a leak opened and the pressure abruptly rose to about 5 p of 
Hg. Since all attempts at achieving a high vacuum in the model had failed it 
became necessary to run the tests under conditions where the gas conduction is 
significant . This meant that the effect of gas pressure on the thermal response 
of the model had to be determined. The pressure must be high enough for the 
gas heat conduction mechanism to be described by the thermal conductivity of the 
gas (as. opposed to molecular or slip flow considerations) ; however, the pressure 
must also be low enough to eliminate any free convection effects. Also the 
increased pressure in the model necessitated expansion of the thermal math model 
to account for the gaseous condition process. 

II .5.6. 1.2 Initial Test Series 

The initial test series consisted of 17 runs including five runs made during 
the "bake out" of the system. These runs during the system "bake out" were 
primarily made to monitor the model temperatures and to checkout the repeat- 
ability of the measurements . The remaining runs were made to determine the 
pressure effects on the model temperatures and to generate the data required for 
upgrading the thermal math model. Table IT-7 gives the test conditions for these 
runs . 

The power input was determined by measuring the voltage and current (by means of 
a calibrated shunt) with a null volt meter. The D.C. power supply regulated 
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the voltage to within about +0.01 volts. The pressure was controlled by 
throttling the flow to the vacuum system and bleeding in air between the model 
and the throttling valve. Generally, several sets of data were taken to 
establish that steady state conditions had been achieved, however only one run 
number was recorded. The stability and repeatability of the experimental con- 
ditions proved to be excellent. After taking a set of data the thermocouple 
output that was measured first was remeasured and if it had changed more than 
1 p volt ( 0.03°F) the test run was continued and another set of data taken at 

a later time. The following table gives the maximum temperature change and time 
elapsed between the last two sets of data taken for the runs. 

Maximum Temperature Change 


Run No. 

Elapsed Time (hours) 

yVolts 

1-6 

Not available 

8.0 

7 

3.0 

8.0 

8 

19.0 

7.0 

9 

7.0 

42.0 

10 

6.0 

5.0 

11-13 

Not available 


14 

5.0 

5.0 

15 

4.0 

15.0 

16 

4.5 

7.0 

17 

3.5 

4.0 


Runs 6-14 were made to determine the effect of gas pressure on the thermal 
response of the model. Runs 6-10 were made at one heating rate and 11-14 at 
a different heating rate. Run 6 was made at the lowest pressure that could 
be attained in the model. After taking the data for this run, the pressure in 
both the inner and outer cylinders was slowly increased while the temperature 
at the heater section (TC 6) was monitored. It was expected that as the 
pressure increased and the heat transfer mechanism for the gas changed from the 
molecular or slip flow process to the continuum process, the heater section 
temperature would decrease because of the increased conductive losses. Figure 
11-37 shows the unexpected results of this test. Instead of decreasing, the 
temperature increased for about 20 minutes before it began to drop. This 
indicates that the thermal conductance between the heater core and the cylinder 
wall is a strong function of pressure (as would be the case if a gap existed 
between the heater and the wall). At the lower pressure the thermal conductance 
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Figure 1 1-37: TRANSIENT RESPONSE OF HEA TER SECTION 
DUE TO CHA NGING GAS PRESSURE 
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is small and the temperature of the heater core is considerably higher than 

v 

the equilibrium wall temperature and as the pressure (and thermal conductance) 
is increased the excess thermal energy stored in the heater core is transferred 
to the wall giving an increase in wall temperature while this excess heat is 
being dissipated. 

Figure 11-38 shows the temperature distribution for runs 6-10. These data show 
the effect of gas pressure on the thermal response of the model. The following 
observations can be made of the data. 

o Free convection effects are readily apparent in run 9 where the outer 
cylinder was at atmospheric pressure. 

o The data for runs 7, 8 and 10 group together with the exception of TCs 7 
and 19. 

o The low pressure case (run 6) shows a significant difference from runs 
7, 8 and 10. 

Figure 11-39 shows the temperature differences between runs 6, 7, and 8 and 
run 10. The temperature for run 6 is significantly higher for the inner 
cylinder and is somewhat lower for the outer cylinder. This indicates that there 
is less heat transfer by gas conduction and the pressure (5 p of Hg.) is low 
enough for the gas conduction to take place by molecular or slip flow phenomena. 
The temperature for runs 7 and 8 are slightly higher at the heater section and 
slightly lower for the rest of the model than those for run 10. The difference 
in temperatures for the inner cylinder can be explained by the higher gas 
pressure in the inner cylinder for run 10. This higher pressure gives increased 
thermal conductance through the microquartz insulation in the inner cylinder 
and lowers the heater section temperature while increasing the temperature in 
the surrounding regions. The generally lower temperatures for runs 7 and 8 
might be explained by a combination of the increased heat leak from the heater 
section through the heater leads, the very slightly larger heat input for run 
10, the differences in the ambient temperature and the differences in cooling 
water temperature. For purposes of upgrading the radiation-conduction aspects 
of the thermal math model, the temperature gradients are of greater importance 
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Figure 1 1-39: TEMPERA TURE DIFFERENCES BETWEEN RUNS 6,7 & 8 VS RUN 10 
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than the absolute values of the temperature. Table I 1-8 gives a comparison 
of the temperature gradients for runs 6, 7, 8 and 10, With the exception of 
the gradients that involve TCs 7 and 19, the agreement is quite good. It should 
be noted that the thermocouple outputs were measured to the nearest microvolt, 
except for TCs 37-44 where the temperature gradients are small. The outputs 
for these TCs were measured and read to the nearest 0.1 microvolts with an 
approximate accuracy of +0.1 microvolts. 

The results of runs 6-10 show that a pressure in the outer cylinder region on 
the order of 1 mm of Eg. will result in gaseous heat transfer by thermal 
conduction. The results also show the temperature distribution for the inner 
cylinder to be dependent to some extent on the pressure inside the inner 
cylinder. 

Runs 11-14 were made at an increased heating rate over runs 6-10 and represent 
a further investigation of the pressure effects. Run 11 was made at a pressure 
of 1 mm of Hg. in both inner and outer cylinders and run 14 was made at a 
pressure of 1 mm of Hg. in the outer cylinder and at atmospheric pressure in 
the inner cylinder. The temperature comparision between runs 11 and 14 is 
very similar to that between runs 7 and 8 and run 10. Run 12 was made with 
both inner and outer cylinders at atmospheric pressure. The free convection 
effects were easily seen in this run. Run 13 was an attempt to run with the 
inner cylinder evacuated to about 1 mm of Hg. while the outer cylinder was at 
atmospheric pressure, however, the small pumping line to the inner cylinder 
and a leak between the inner and outer cylinders limited the pressure to about 
3-5 mm of Hg. Aside from the pump down problems, refrigerator problems during 
this run gave rise to variations in the inlet cooling water temperature of about 
1°F which invalidated the data taken. 

Since it was apparent that a reduced pressure could not be maintained in the 
inner cylinder for the cases to be run with the outer cylinder pressurized, it 
was decided to run the radiation-conduction test series with atmospheric pressure 
in the inner cylinder while maintaining a pressure about 1 mm of Hg. in the outer 
cylinder. Runs 10 and 14 constitute part of this test series and runs 15-17 
make up the rest of the test series. 
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TABLE I 1-8 

TEMPERATURE GRADIENT COMPARISON 


THERMOCOUPLES 


TEMPERATURE GRADIENT 

* 

(y volts) 



Run 6 

Run 7 

Run 8 

Run 10 

1-2 

391. 

385. 

388 . 

393. 

13-14 

TC # 14 NOT 

READING 



2-3 

841. 

782. 

784 . 

785. 

14-15 

TC # 14 NOT 

READING 



3-4 

1057. 

1012. 

1013. 

1010. 

15-16 

1041. 

1010. 

1011. 

1007. 

4-5 

1445. 

1376 . 

1373. 

1359. 

16-17 

1415. 

1381. 

1384. 

1371. 

5-6 

1140 

1097 

1098 . 

1083. 

17-18 

TC # 18 NOT 

READING 



6-7 

1025 

516 

447. 

329, 

18-19 

TC # 18 NOT 

READING 



7-8 

2197. 

1631. 

1563. 

1420. 

19-20 

1524. 

1312. 

1290. 

1230. 

8-9 

1480. 

1401. 

1399. 

1388. 

20-21 

1455. 

1416 . 

1420. 

1417. 

9-10 

1072. 

1026. 

1027. 

1022. 

21-22 

1044. 

1022 

1025 

1027. 

10-11 

840. 

802. 

801. 

801. 

22-23 

834. 

805, 

808. 

809. 

11-12 

TC # 12 READS ERRONEOUSLY 



23-24 

373. 

374 

376. 

381. 

25-26 

143. 

135 . 

135. 

138. 

26-27 

228. 

220 . 

220. 

223. 

27-28 

228. 

220. 

220. 

223. 

28-29 

323. 

317. 

317. 

319. 

29-30 

443. 

441. 

442, 

445. 

30-31 

5. 

8. 

7. 

9. 

31-32 

425. 

427. 

428. 

431. 

* 35 pv 55S 

1°F 
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TABLE 11-8 (CONT.) 


THERMOCOUPLES 


TEMPERATURE 

GRADIENT (p volts) 



Run 6 

Run 7 

Run 8 

Run 10 

32-33 

321. 

318. 

319. 

320. 

33-34 

226. 

222, 

221. 

224. 

34-35 

168. 

164. 

165. 

165. 

35-36 

120. 

128. 

129. 

130. 

37-38 

46. 

45.6 

46.0 

46.5 

38-39 

34. 

34.6 

34.5 

34.3 

39-40 

28. 

27.2 

26.9 

27.2 

41-42 

64. 

59,5 

58.5 

59.1 

42-43 

45. 

41.2 

41.0 

41.1 

43-44 

32. 

30.4 

30.6 

30.8 

45-46,47 

778. 

787. 

790. 

794. 


386 



D 180- 150 48-1 


As a point of interest, the systematic variation of the temperature around 
the inner cylinder is shown in Figure 11-40 for various pressures and in 
Figure 11-41 for various power inputs* The difference in temperature from one 
side of the inner cylinder to the other is larger than expected. An adequate 
explanation of these differences has not been developed, however, they appear 
to be connected to the heater installation and the heater lead routing. The 
most perplexing phenomena is the reversal of the temperature differences for 
the low pressure cases (run 6 in the data presented, runs 1 and 2 also show 
this reversal). The effect of power input is much greater for the end of the 
inner cylinder (TCs 1-6 and 13-18) through which the heater leads are routed than 
for the other end. The effect of these anomalies will be minimized by using an 
average temperature in the upgrading of the thermal math model. 


II. 5. 6. 1.3 Free Convection Test Series 


After completing the initial test series the 1/4 scale model was set up for the 
free convection test series. The test setup for these runs is shown in Figure 
11-42. The test setup was changed between the 1/2 and 1 atmosphere pressure 
runs and the 2, 4 and 8 atmosphere pressure runs. A photograph of the free 
convection test set up is shown in Figure 11-43 . 

Some difficulty was encountered in running the 1/2 atmosphere pressure cases. 

The pressure was set by adjusting the bleed and throttling valves shown in 
Figure 11-42. The pressure was very sensitive to small adjustments of these 
values. After setting the valves to obtain the proper pressure the system 
was allowed to come to thermal equilibrium, however; during the time between 
setting the valves and taking the data the pressure would change slightly. 

During the first attempt to run at 8 atmospheres an air leak between the inner 
and outer cylinder was discovered. The 8 atmosphere pressure probably increased 
the size of a small leak between the two regions. In order to continue the tests 
of the model a bracket was constructed to retain the air distribution plugs in 
the inner cylinder and the inner cylinder was pressurized to the same pressure 
as the outer cylinder. This support bracket may be seen in Figure 11-43 . After 
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Figure 1142: 1/4 SCALE MODEL FREE CONVECTION TEST SETUP SCHEMATIC 
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making the first run at 8 atmospheres pressure (run 30) a repeat run was made 
for a 4 atmosphere pressure case. Figure 11-44 shows the differences between 
runs 29 and 29.1, without and with the inner cylinder pressurized, respectively. 
This figure indicates that the air leak was removing some heat from the system 
since run 29 has a lower average model temperature than run 29.1. The 
depression in temperature in the heater section region of the model for run 29.1 
is due primarily to the increased free convection (run 29.1 had a pressure of 
4.22 atm whereas run 29 had 4.01 atm pressure). All of the 4 atmosphere pres- 
sures cases were rerun with the inner cylinder pressurized. A check of the 2 
atmosphere pressure cases (runs 31.3 and 37.2) showed no significant difference 
between inner cylinder being pressurized and not pressurized. The free convec- 
tion test series conditions are given in Table II-9. 

II.5.6.1.4 Forced Convection Test Series 

The Forced Convection Test Series test setup for the 1/4 scale model is shown 

in Figure 11-45. Figure 11-46 shows a photograph of the test setup. Calibrated 

"rotometers" were used to measure the air flow rate through the model. The air 

pressure at the rotometer was regulated to a value P q . The system pressure 

P and the air flow rate were set by adjusting the "flox* rate control valve” 
s 

and the "system pressure control valve.” A mercury manometer was used to measure 
the pressure drop AP through the model. This pressure drop was significant for 
the high flow rate cases. Thermocouples were installed in the air stream at the 
ends of the air inlet and air outlet tubes. For the 1/2 atmosphere pressure 
runs the vent line was attached to a vacuum system (the ”roughing” pump for 
Boeing’s Space Simulation Chamber A). The vent line and "system pressure control 
valve" were removed for the 1 atmosphere pressure runs. The operating conditions 
for the 1/4 Scale Model Forced Convection Test Services are given in Table 11-10. 

The thermocouple output stability obtained in the "Radiation-Conduction" and 
"Free Convection Test Series" was not matched in the "Forced Convection Test 
Series.” Temperature oscillations of up to about +0.5°F were present for some 
regions of the inner cylinder and for the air outlet temperature. The air inlet 
temperature variation during a recording of the test data was generally less than 
about +0. 1°F. The outer cylinder and end plate temperatures were quite steady 
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Figure 1144: TEMPERA TURE DIFFERENCE WITH INNER CYLINDER PRESSURIZED 
(DIFFERENCE BETWEEN RUN 29. 1 AND RUN 29) 


393 



FREE CONVECTION TEST SERIES 


D 180 - 15048-1 


Cfl 

CO 

00 

vO 

1 — 1 

co 

CM 

in 

in 

sfr 

rH 

1—1 

o 

M3 

00 

cp 

00 

m 

rH 

oo 

rH 

vD 

CP 

CM 

CO 

i — 1 

CO 

CT\ 


U 

1 — I 

CM 

p^ 

-<f 

p. 


cp 

CP 

rH 

o 

r-i 

rH 

CP 

<p 

vD 

H- 

Mf 

p. 

<p 

CO 

CO 

CM 

O 

O 

p* 


co 

mT 

4-» 

00 

00 

r- 


• 

• 

* 

.* 

« 

• 

• 

• 

• 

• 

• 

p~ 

p» 

• 

* 

• 

• 

-• 

• 

* 

• 

» 



cd 

• 

• 

• 

• 


< — i 

CO 

co 

GO 

00 

00 

00 

CO 

CO 

i — 1 

.• 

• 

t — t 

co 

00 

00 

oo 

<r 


rH 

rH 


* 

56 

m 

m 

00 

00 

H 

rH 

T — 1 

r — 1 

rH 

iH 

I — \ 

rH 

rH 

iH 

iH 

00 

00 

rH 

i — 1 

»H 

rH 

1 — l 

«H 

rH 

iH 

rH 

00 

00 



4J 

a 

< 

cp 

CM 

<N 

p* 

p- 

m 


p« 

CM 

p- 

o 

o 

CP 

CM 

CM 

o 

O 

vO 

CO 

CO 

o 

CM 

a\ 

O 

00 

o 

00 

o 

w 

a) 

X 

rH 

CM 

O 

CP 

o 

o 

p> 

p. 

CO 

o 

CO 

CO 

p- 

00 

O 

o 

o 

O 

00 

in 

MO 

m 

00 

Ov 

o 

rH 

CP 

o 

J3E 

u 


rH 

rH 

vO 

m 

o 

o 

CM 

CM 

P» 

p» 

p* 

p. 

CM 

CM 

o 

vO 

vD 

O 

CM 

ip 



CM 

CM 

o 

o 

m 

vD 

o 

u 

< — i 

CM 

CM 

CM 

CM 

CO 

co 

CO 

co 

CO 

CO 

co 

CO 

CO 

CO 

CO 

CM 

CM 

CO 

CO 

CO 

CO 

co 

CO 

CO 

CO 

CO 

CM 

CM 

0-1 

d 

o 

• 

• 

* 

• 

• 

• 

• 

• 


* 

* 

> 

• 

• 

.* 

• 

• 

• 


• 

• 


• 

# 

• 


• 

• 


a 

o 























Q 

o 

o 

o 

o 

o 


a a 

o o 

•H *H 

4-3 LO 

vo cp n3 vo vo oo h h rtin^r^^H^rrNOONoofocn^ooooooooNrs 

<t'd‘*HP^vDOOvOvO*HvOvOminvOvOavvOvOCrivOOOOOaOvOvOavOvOvO 

* * j_i • ••« • • M • * • • • • • • • • • • • « « • « * • • 

P* P* ctf CO CO Ov Ov CM CM CdOOOOOOOOCNCMOOCOcOOOCMcOOOOOCMCMOOavOOcO 

cm cn > co co co co <r ^ st <r <r mj- ro co co co st sr >d- <r -<r co co co co 

co m <r vo h m oo co cs .in 

HHHrlOOiOHON <fOOO\v£)vO\0(OOlHOCr»OHCNH^OcnHcO 

m o o o <r sr o o sr ^<j\.^b.o0ooffioaiOpNo.oH^oooi 

OHHHOOHH.O O O O M d N N CS (O ^ CO ST 00 ^ Sf ^ N <f <T N 


sr m oo < j \ lo o 

«— i ^ Mf P> -d vO vO m O H 00 CM CM CM • CM H • • vD CM • 

• •• • • • • *••••«♦• CM •• rH rH • • CM 

in m \o in to o n in :.n n m cs ^ sr o ^ ^ o o ^ o 

rHOOOHiHOOrHr- 1 rH o o i l *H I— 1 rH -d* <t -d- ST rH <1- -sf rl ri «d" *d" i — I 


tnr^<tovcooNvoco 
<r CO CO co O M H N 


•d-pscomo>d-<rp^vDLnvDrHP>rHOO^ointHavp^ 

vOcooooooovTHvoc^r^.cooO'd-cocNcoavr^.aNr^ 


O O O O O ov O O ov 

CO CO CO CO CO M CO co <N 


ovaNOvovoovoovc^avovovooooovaNOvov 

NMNOJcONCOCMNiNCNCNCOcOfOcOcNCMCNCM 


a) 

O 

m 

O 

O 

O 

O 

O 

o 

o 

O 

o 

O 

O 

o 

o 

O 

O 

O 

O 

o 

o 

o 

m 

m 

m 

o 

o 

O 

S 

CO 

r — 1 

CO 

CO 

CO 

CO 

o 

o 

O 

o 

CO 

CO 

CO 

o 

CO 

CO 

CO 

o 

CO 

o 

CO 

CO 


rH 

<3* 

o 

CO 

O 

•H 

av 

in 

P- 


rH 

<r 

rH 

00 

00 

00 

<* 

vO 

O 

m 

rH 

00 

o 

V0 

00 

00 

m 

m 


vD i 

1 <t 

00 

<r 

o 

H 

o 

rH 

rH 

i — 1 

rH 

rH 

rH 

rH 

»H 

1 — 1 

rH 

iH 

i — 1 

rH 

r — i 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 


CP CT\ CP CJV CPCPCPCPCPCPCPCPCPCPO'vO'vCP 

vO vO vO vD CP CP CP CP Ov CP vO vO vO vO vO vD vO vO vO vO vO vO ^D OO 

" — \ \ \ \0 vO vD vO vO vO *\ '>*- - — » Ss, p** \ \ Ss *>>. \ \ \ O O O O O P^ P-- 

co n w'x's wo h cl n n n oo o> OiHcMcotHf^r^r>r^r^>->.>< 

H^CNCNHcOsfinOOa^HHHHHHHHNCMCMOJCOWWNOCNl 

HHHH(NCMNCNiMCN}CN|CM(NCMNCM(NCslCNCMCNjCN<MW'X\ v vW 


oo^oOHnjco^in^^NNOOHH^MXJinmc^o^cNcovoinoo^ 

HHMCMCMcMCjCMcMNCMCStCNlCMCOcOcOCOcOCOCOMCOCvlcOcOcOcOCOCO 


394 



TABLE II-9 (CONT 


D180-15048-1 


CT\ 


cn 

CO 

• 


co 

* 

i — l 

00 


CM 




00 

CTk 



CM 

m 



co 

CM 



• 

* 



o 

o 


CO 



to 

s 

m 

in 

(U 

d 

■V O 

vD 

a 

u 

mo 

vO 

•H 




cd 

u 

CM 

CO 

+J 

0) 

St 

CO 

d 





0) 


ON O 


<D 


<J\ o 


40 a) 


m • 

CO 

40 


rH CM 

d 

o H 



d 

u 



u 





H • 



u 

a) co 



Q) 

TO d 



40 

d d 



4J 

•H U 



O 

rH 




Po 5P 


uo 

u 

a a 



o 



• .* 

MH 

u £ 


sr <r 


CO o 


i— i i—i 

00 

d i— 1 



•H 

d th 



CO 

•H O 



PM 

MS 




TO 

0 


a 

a> <0 

a) 


•H 

U 40 




•H H 

a> 


TO 

d 

40 


a 

cr to 

CM 

00 «>± 

cd 

<0 d 

CO 

m 


u cd 

o 

• « 

r^. 


6 

cri 

CM 

M o 

jj 

CM CM 

1 

a) co 

cd 


00 

TO 



1 — 1 

d d 

0) 



■H d 

40 


CO 

t— 1 M 

u 


d 




d 

U M 

o 



O 

4J 



H MH 



u 

<0 

d 


o 

4-» <0 

0) 


M-l 

d m 

PM 



o d 

o 


00 

CO 



m 

TO CO 

u 



d <o 

<u 


MS 

cd H 

TO 


o 

PM 

d 

O O 


H 

•H 

CO o 

CO 

CO U 

rH 

O MO 

<u 

d <u 

Ps 

i—l rH 

40 

d TO 

a 


a 

•rl d 



a 

•H 

u 


•H 

d i— i 

a> 



CD Po 

d 


a 

cu o 

d 

o o 

•H 

£ 

•H 



■U M 


\ S* 

00 

0) 0) 

0) 

<r m 

a 

40 4J 

40 

i—i r—i 

•H 

d 

4J 

\ \ 

TO 

4*5 O 


T— 1 »H 

cd 

cd 

TO 


<u 

a) -u 

cd 



Hi d 

40 

CM *H 


* 


• • 

-K 

* 


i — i r> 




CO CO 





395 



SYSTEM 


D180-15048-1 



396 


Figure h 









^4Ui. 




spill; 








..'.... 


tlSSS 


ifpS 

i&i® 


FIGURE 11—46 FORCED CONVECTION TEST SETUP 







TABLE 11-10 1/4 SCALE MODEL FORCED CONVECTION TEST SERIES 


D180-15048-1 



Pd 03 


00 on O NO CO 

co on oo oo 

oo in no 

St ON 

st 

co 

ON o 


rH C3 ON ON NO S* NO 00 

rH CM H rH 

S' ON ON 


M H 


00 cO O ON o 

CM rH CO CO 

co co co 

oo S' 

00 

00 

ON o 


<t cn m ch a\ 

ON ON 

rs S* S' S' 

st st sf 


3 H 

o 












O <3 


oo oo co s. oo 

oo 00 oo oo 

oo oo oo 

S' S' 

S' 

S' 

co st 


COOOOOQCCOCOCOCO 

HrlrHHO 

o o oo oo oo 


Pm 3 


i— 1 rH CM CM CM 

CM CM rH rH 

H H H 

CM CM 

CM 

CM 

rH rH 


H H H HH H 

rH rH 

HH HH 

rH rH rH 

H 

03 

& * 














© 

JS 














PH 

<! 














a 

H M3 














w 

S X 00 03 

















r» n N 0\ M 

rH oo m m m o h m o st 

co 

m no 


M3 st CM rH H CM 

CM St 

00 S' S' NO 

m S' S' 


S H cr> s 

o 

if) nO st 00 St 

NO in NO NO NO NO nO 

CM CM 

CM 

CM 

oo oo 


M3 nO NO NO 00 00 

oo oo © o o © 

m S' S' 


Bd o o <3 


N S vO \0 iO 

NO NO s N 

i — r-. r^- 

NO NO 

NO 

NO 

CH CM 


S S N S N tM 

CM CM 

© o o o 

S' S' S' 

o 

B °-o 


CO CO St Sf St 

St st CO CO 

co co co M -t -t 

st 

co CO 


n ro n n co m 

**•••« 

CO co 

co co co oo o 

O O CO CO co 


rH V* 


O O O O O 

o o o o 

o o o 

o o 

O 

o 

o o 


o © © o o o 

o o 

o o o o 

o o o 


8 H 




St rH 


m 


co no 


no m co S' 

CM 

CO S' st 



<5 © 


in n oo m H 

co oo n. rs 

n> in no 

m cm 

S' 

m 

CM CO 


CO CM CM co o O 

ON © 

ts o S' 00 

00 ON O 


H O 

o 

co oo co co <r 

no m oo oo 

oo oo oo 

CM CM 

CM 

CM 

NO S3 


ON ON ON ON NO NO 

m M3 

ON o ON ON 



© > 















o 


oo co o O O 

o o po oo oo oo oo 

o o 

o 

o 

CM CM 


00 00 00 00 CM CM 

CM CM 

00 on 00 00 o 

o o oo oo on 


> 


Sf sf N0 NO NO 

no no m sr m -t 

NO NO 

NO 

NO 

sf St 


•fstst'tsfStstst 

co co co cl 

st St st 


a 


ON On ON H N 

in 

00 




NO 


Ss 






S' 

m so s. S' 

S' S' o 

O r>> 

m 

o 

CM f> 


ON NO oo On h on cm m 

CM CM H M3 St 

m co cm m on 


W 03 

ON 

ON On O' H 00 

00 ON ON O 

o ON o 

O O ON 

© 

o o 


03 O On ON O On O © 

o o o o o 

ON rH rH O ON 


H 03 


• « * • • 












03 H g 

rH 

rH rH rH rH CO 

co s S H 

HHM00 H 

t— 1 

st 

St rH 


S'SfrHOrHrHstOO 

00 rH CM St CM 

CO rH rH CM CO 


S H 

03 <! 















o 















a 

















in in m m m 

in m m m 

mom 

m m o 

o 

o o 


m 


m m 

m m m in in 


ru ph 

00 

S' S' s- o st 

sf NO O CO 

C0 CM H 

O CO CM 

rH 

rH CO 


o 


© S' St CM CM S' st 


< o 
















CM 

CM CM CM m rH 

rH O O O 

o o o 

o o 

O 

O 

o o 


o 


O CM 

rH m m CM rH 


a 









Vi 



VH 











Vi 


















o 














a 

a 



in 

m 

m 



m 






a 

03 

03 


Pm 

030 O 

tu M 

O 

co co co on in 

cm in m cm 

CM CO CM 

O CM 

CM 

CM 

m cm 


in CM S' CM CM NO 

co m o no oo .oo 

m O O rH CM 


m 

m m m cm cm 

CM CM CM A 

H M M 

CO rH 

st 

st 

St rH 


tHmst©©st«tcocoostsrm 


a 


03 2 

rH 


-too 

rH St 

o 

rH 

St 

st 


O St rH rH St o 


st rH St 


p-< a 



rH rH 


rH 





rH 



ph 
















O ,-N 















M O 

Pd o a 















h a 















§ r 















o o ° 

ON 

<{ .co .<! < <sj 

cm O <J IS. 

ON O O ON <ri 

O 

NO 

S' S' 


ON S' M3 00 S' H 

CM rH 

CM in St CM rH 

rH O S' CM CM 


pd _ a 

m 

-s. U0 

no no "Sm 

m no no 

m m 

m 

m m 


St st St st m NO 

NO S' 

S' S' S' S' S' 

S' S' nO NO nO 


< a m 

S' 

a n* 52 a a 

s« s. a S' s> s. s. 

S' a 

S' 

S' 

S' S' 


S' s* S' rs. S' 

S' S' 

S' S' S' S' S' 



« S'-' 














a 















o 
















a 










St CO CO ON O st 
co st <t st m st 

ON O 

m CM CM rH 






o CM 

m s- cm 

CM St 1 

m 

CM 

CM CO 


co m 

in m m no 



a 

in 


NO NO ON ON 

co oo on 

ON ON 

00 

ON 

ON ON 


S' S' S' S' S' f — 

S' S' 

S' S' S' S' st 

st st st st sf 


a 

CM 

CM CM CM CM CM 

CM CM CM CM 

CM CM CM 

CM CM 

CM 

CM 

CM CM 


o o o o o o 

o o 

O O O O CM 

CM CM CM CM CM 


3 as 

rH 

rH rH rH rH rH 

rHrHOOOOOOO 

o 

o o o 


0000©©0©0©00rH 

rH rH rH rH rH 



o 

a 










mcMmomrHom 

o co o m 




H 

m 

o CM 00 00 o 

NO CM O O 

S OO o 

o © 

S' 

o 

o o 


00 © O CM CM rH 

O CM 

cococostooocMcMinm | 




M3 

no no in m no 

ON ON CM CM 

H H N 

CM CM 

rH 

CM 

CM CM 


00 ON ON ON ON ON 

ON ON 



w 

H 


a 

m 

m uo m m in 

CO CO rH rH 

rH rH rH 

rH rH 

t-H 

rH 

rH rH 


rH rH rH rH H H 

rH rH 

rHrHrHrHmmmmmm 

a 

a 

H 

o 














o 

g 

t— i 
03 
P-t 














£ 

O 

o 

o m m cm oo 

o cm in o 

o o o 

o o m 

O 

UO O 


m co o m m cm 

m m 

oo©moooooooomm 

Pd 

H 
















Q 

* 

rH 

o o o o o 

ON ON 00 o 

o o o 

O CM 

O 

o 

O rH 


on on on oo oo oo 

oo oo 

OOOOOOCOOOrHHOO 

M 

Pd 

o 

m 

m m m m in 

O O O rH 

rl H H 

rH rH 

rH 

rH 

rH rH 


o © © o o o 

o o 

o o o o m 

m in m m m 



PP 



rH rH >H rH 

rH rH H 

rH rH 


•H 

rH rH 

CO 

H H rl rH H rl 

rH rH 

rH rH rH rH 



' . 

r- 


CM 






~T~ 

CM 



S' 




o 

M3 


O 






CM 



NO 




a 

ON 


rH 






m 



ON 





in 


NO 






oo 



m 





o 















g 

m 


H 






rH 



m 




Pd 

q 


X 






q 



o 




m 

< 


u 






<3 



< 




03 

m 


« 






PQ 



Pd 




o 

o o o o o 

o o o m © © © 

moo 

o o 

O O O 

©mo©©inmo©ooo 

m o o o o 



o 

CO o o o o 

o o o iH 

© CO CO sf © 

m 

O CO 

CO CO 

CO 

CO rH CO O CO 

32 

coocoocostcooom 



oo 

o in oo co r-n 

NO co rH Sf 00 CO S' 

co on 

NO 

rH CO 

CO H O' 

st ON CO ON st 

st CO oo co o 

rH CO O m S' 


EH 

rH 

rH rH rH CM rH 

rH CM CM rH 

rH CM H CM O 

rH 

rH rH 

rH CM 

Q 

hhmohh© 

rH CM O rH CM 

CM CM H rH rH 





o 

O 

o 


© 


o 

o 

o 

o 

© 



O 

O O 

o s. 

f> 

S' 


S' 


S' 

S' 

S' 

S' 

S' 



IS. 

S' fS 

s. \ 







'S' 



S' 


W 



"s. o 

H 

CM 


CO 


NO 

S' 

oo 

ON 

o 


Eh 

St 

in m3 

ON rH 

rH 

rH 


rH 


rH 

rH 


rH 

CM 


<3 


S'. >s 



S'. 






S' 




Q 

CM 

CM CM 

CM CM 

CM 

CM 


CM 


CM 

CM 

CM 

CM 

CM 




rH 

rH 




rH 








53 • 















5 o 

CM 

n co sr m no 

NO S' 00 ON 

O rH CM 

CO St 

m 

NO NO 

s* oo 

ON 

© rH CM CO St 

m \o 

S' CO ON O rH 

cm co st m no 


a a 

st 


t t MM 

m m n 

m uo 

m 

m m 

m n m 

NO NO NO NO NO 

NO NO 

M3 NO M3 S' S' 

S' S' S' S' S' 



ai 

cd 

© 




o 








398 




TABLE 11-10 1/4 SCALE MODEL FORCED CONVECTION TEST SERIES (continued) 


D180- 15048-1 


1 a 

Pu ? 


1 — I tH tH rH P- .ON VO ■ St 
ts-ococoonvocoon 


PH 

S5 



mmstOONrHvOvO 
OWiniflNOinN 
O cm rs. ts. cm © rs. cm 
eocococococococo 


oooooooo 


w 

q 

<J cn 

E-t H 

O O 
> > 


H H cm n co h in 
OOOtHOOOO 
in in m in in m w> in 


oooooooo 


Pd 

pm 

< PH 
O 


<1 

a 


o 

a 


O Pm 

U)H O 


HHcnOiCM\OH-tf 

OOOOCOOOCOOOOOOO 

cococncncocnnn 


u ^ 

M o 

B 

1 

Pd 2 


<J <5 <3 

• — o o cm oo ro 

Z Z 2 ^ CO CM o ON 


8 


a pq 

t-j 


CONHM 

IS in in 00 N N H M 

lO MO O St St M 00 00 

|s. t-. r>. o © © »H rH 

OOOOOOOO 


OOOOOOOO 


O 

a 




LHCOOOOOOO 


o 

:S 


o 

a 

a 

s 


00 _ 

St 

tH 


1 


oomoooom 

Cl rl St CM CO O CO St 
cors.oois.cM on is«o 
OrHOrHCMOrHO 


m 

1 


o 

is. 

sr 


CM 

cm 


o 

Ss* 

M3 


O 

a 


NOOOiOHCMCOst 
N S is 00 00 CO 00 OO 


H H 00 00 00 00 oo 


o o 


oo -t o M) n m n 

<t- st m st -3- st st 

oo oo co oo oo oo oo 

rH rH rH rH rH rH *H 


CM CM ID o in 

CO St |S fs \£> (ft OO 
O O f-s rs ts r-. 
O O co cn co cn co 

CO CO o o o o o 


m oo o in fs H co 

Js. l£> CO fs, !>• fs* fs. 

p- r — is p> is. in ts 

co co co co co co co 

© o © © © © o 


h m p- rs st on m tH «h 

oiooist ocM in o\ vd oo co co m co 

Oi O ON cti ov ct> OI OO 00 00 CO 00 OO OO 

OOOVCOOOOOOOOOOOOOOOOOOOOOOOOO 
co CO -St St sfr -4- st st st st st st St st 


cn is 00 CM VO ON -00 VO 

©ONCOCMCMHincMsfCOCOHstONOON 
instHOOiOONOHOOOOOMinst 

OOH'tioHNHHtvtHCMOOO 


<3 

"s. in m is. in 

a PCI W CO o cm fs. o 

0 0 0 0 0-3-0 

V 1 Vi 


3 ss a a 

o 

jlL 


/-N CS /s /s 

cm co o\in n cn o (O cm o n o on mj is 

00 IS • fs. oo IS. 

cococMstcOocMocorHstomcococo 
s_^w st St o 0-3- Hs/s/v 

rH .rH 


<1 /-N / — \ /-s /-s / . /-N /S /*S /S rs <J 

Sv St MO St CM CM mmsf lOOCM CMNiOS 

aPOvOvDvOvOvOvOvDvO is is. |s J2J cM CO 
• is.rs.r^r>r^r^.i>.r^rs.|s.r^ • • 

OV W v — % — — / s-' S^ O O 

CM CO CO 


rH ON o O -3" st St CO 00 

oors. co st st co st m st st st co st 

HHststONONONONOOMsMsrsHO 
OOCMCMCMCMCMCMCMOOOOOOO 


o O rH rH o O © O *H O O O © O © © 


mm mm 

onOOOOOOOOOONOOOOOOO 

n in mo vo cm cm cMcMstcoovONONONvom 
HHmmHrlrtHmHHrlHHH 


OomooooocoomooooorH 

OOOrHOOONONOONrHOOOOO 
m m i — I rH O O rH O rH rH rH rH 
iH tH rH tH H iH iH iH rH rH 



CO 

fs. CM 

CM 

vO O 

CM 

ON rH 

m 

m so 

pr 00 . 

O fs. 

St 

m rH 

rH 

Q X 

U 

< u 

< 

pq PQ 

PQ 


© © o © © © m o © m m o m © © o 
OrHOCnocOstOOtHstCOstcOCOO 
oiinstcooMnoNstocointoocMON 
ohhhhhcmohoohnhcmo 



CO CO CO cO CO CO 


iH 

mmvors.oooNO H cm rostmvo^oooN 

00 00 00 00 00 OO ON ON CJN CJN ON ON Os ON ON ON 

no 




tj 



M 



cd 

1 



I 

T3 


u 



o 

> 


o 

O 


TJ 

rH 



IW 


a 



d 



tH 

a) 

DO 

rH 

CO 

O 

u 

PI 

u 

CO 

•H 


J> 

4J 

00 

0) 

4-1 

d 

Is 

(0 

*H 


CO 

T3 

d 


rH 

o 

u 

•H 

*H 

(0 

d 

4-1 

4J 

a 

O 

si 


(0 

a 

iw 

Is 

o 

o 

•H 

4J 



o 

00 


Is 

d 

> 


*H 

o 

4J 

d 

tH 

ti 

a) 

VH 

(0 

cu 


u 

o 

H 

<U 


•H 

4S 


d 


*Q 


•H 


a) 

a a 

4J 

> 


0) 

cd 

Xi * H 

CO 

a 

H l-l 

CU 


•H ,D 

d 

ON 


S 

ON 

rH 

1 

d *h 

d 

vO 

d d 

«H 

oo 

u cr 

1 u 


cu 

X 

CO 

4-1 

•H 

d 

Cfl 4-1 

rH 

d 

<0 cd 

*rl 

a 

a 4-i 

d 


<0 o 

cr 


a a 

w 

cu 

4-1 

* • 

• 

O 

rQ O 

T3 

a 


399 




D 180- 150 48-1 


with no measurable temperature oscillations present. At the lower heating 
rate cases the temperature oscillations were much less than for the higher 
heating rate cases. 

11,5.6.2 Full Scale Model Tests 

Prior to testing the full scale model, the thermocouples were checked for 
continuity. Several broken leads were discovered. These were all repaired with 
the exception of thermocouple number 13. Thermocouple numbers 48 and 49 were 
installed in the air distribution plugs for measurement of inlet and outlet 
gas temperatures and the model was insulated with polyurethane foam (as shown 
in Figure 11-31). 

The 1/4 scale model test setup, without the pressure control equipment, was 
used for the full scale model tests. A photograph of the full scale model 
and the test setup is shown in Figure 11-47. 

The "Free Convection Test Series" operating conditions are given in Table 11-11. 
These tests were made with the bottom air supply tube blocked off to prevent gas 
flow through the model. The top tube was left partially open to keep the model 
at atmospheric pressure. 

The "Forced Convection Test Series" operating conditions are given in Table 
11-12. It is interesting to note that the small temperature oscillations 
present in the 1/4 scale model forced convection tests were also apparent in 
the full scale model tests. 

II. 5. 7 Data Analysis and Correlation 

II. 5. 7.1 Thermal Math Model Expansion and Upgrading 

The basic thermal math model for the 1/4 scale model was expanded to correspond 
to the actual test configuration of the model. The nodal model expansion 
included the following: 
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o Foam insulation (polyurethane foam) 
o Insulation inside inner cylinder (microquartz) 

o Support rod (phenolic) 

o Air distribution plug (phenolic) 

o Heater leads (copper) 

o Air flow tubes (stainless steel) 

o Gas inside of model (air) 

Figure 11-48 shows the nodal network for the expanded thermal math model. The 
inner cylinder assembly is represented by 26 nodes exclusive of the inner 
cylinder wall nodes of the basic math model. The air inside the model is 
represented by 28 nodes and the foam insulation by 49 nodes. The external 
temperature (ambient) is represented by a single boundary node. The basic 
model structure is represented by 12 nodes for the inner cylinder, 8 nodes for 
the end plates, 12 nodes for the outer cylinder and 2 nodes for the cooling 
fin (the basic math model consisted of these nodes) . 

The experimental data from test runs 10, 14, 15, 16 and 17 were used in con- 
junction with the expanded thermal math model to determine the upgrading 
necessary to develop a valid thermal math model. The temperatures for these runs 
are presented in Appendix II-A.1.1. These temperatures were calculated from the 
experimental data using the thermocouple calibration correction. The tempera- 
tures for the inner cylinder were taken as the average of the two thermocouple 
(where available) measurements with the exception of node 307. Since the 
thermocouples (numbers 7 and 19) used to measure the temperature at node 307 
apparently are reading high and since the radiation-conduction tests give a 
symmetrical temperature distribution, a temperature correction term was developed 
for thermocouple number 19. This term corrects the temperature measured with 
thermocouple 19 to that measured by thermocouple 6 for the radiation-conduction 
tests. Figure 11-49 shows the correction term as a function of heater power. 

This correction was applied to all of the 1/4 scale model test results listed in 
Appendix II-A. The temperature distributions for test runs 10, 14, 15, 16, and 
17 are shown in Figure 11-50. 
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Figure U-49: TEMPERA PURE CORRECTION FOR TC NO. 19 AT NODE 307 
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Two additional 1/4 scale model radiation-conduction test runs were made (runs 
40 and 41, listed in Appendix II-A.1.1) without water cooling in an attempt 
to obtain a better definition of the model heat leak paths- The results of 
these runs were used in most of the thermal math model upgrading runs. However , 
these results were of limited value and the upgrading effort concentrated on the 
results from test runs 10, 14, 15, 16 and 17. 

The required upgrading was determined by comparing the basic math model thermal 
conductances calculated from the geometry and thermal conductivity (input con- 
ductances) with those calculated from the experimental data (upgraded conduct- 
ances) . The method of obtaining the upgraded conductances and comparing them to 
the input conductances is as follows: 


1. The expanded thermal math model nodal network is input to the BETA program. 

2. A subroutine takes the thermocouple output data from the test run, applies 
the themo couple calibration correction and determines the temperature at 
each node of the basic math model and the ambient temperature. These 
temperatures are then made boundary temperatures for the expanded thermal 
math model. 

3. The BETA program then calculates the nodal temperatures for the expanded 
thermal math model using the boundary temperatures obtained in the sub- 
routine and the measured power input to the heater. 

4. A subroutine than takes the calculated temperatures and makes the following 
determinations for each node of the basic math model. 

a. Radiative heat transfer to the node. 

b. Conductive heat transfer from the air inside the model to the node. 

c. Conductive heat transfer to the node from one (two for two nodes) of 
the following sources, the heater, the insulation inside the inner 
cylinder, the air distribution plugs and the foam insulation. 
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5. The subroutine then makes a thermal balance at each node of the basic 

math model and determines the thermal conductances between adjacent nodes. 
The method of this determination can be explained by referring to Figure 
11-51 which shows the nodal network detail near the heater. The power 
input is divided equally between nodes 526 and 527 which represent the 
heater. The conductors connected to these nodes represent the various heat 
flow paths from the heater. The paths of interest for the basic math model 
are those between nodes 526 and 306 arid nodes 52 7 and 307. By assuming 
negligible thermal conduction between nodes 306 and 307 the thermal con- 
ductances may be determined stepwise around each half of the model start- 
ing at nodes 306 and 307. The radiative heat transfer to node 306 is 
given by 


**>*<**> - f [t*- 


the heat transfer to node 306 from the air is given by 


QGAS (306) - K(306) [ T(606) - T(306) ] 
and the heat transfer from the heater is given by 


QINS(306) = K( 106) [ T(526) - T(306) ] 


The heat conducted between node 306 and node 305 is then given by 
(QRAD(306) + QGAS (306) + QINS(306)) and the upgraded conductance 
between nodes 306 and 305 is given by 


K'(5^ = 


<S?fZbO CSoC,*) + Q<3A SCSOC,') <QIt4S ( 3Cu 


The conductance between node 305 and node 304 can now be calculated 
from 



Q gfrDfgdgj +Q 6 A 1 S £3bs) +OrMSC3oS> Kto D 
T( So*) - -CC3o4) 
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Figure 11-51: THERMAL MA TH MODEL NODAL NETWORK DETAIL NEAR REA TER 
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where 


QRbDC3o&- IS V A -jog T*(SO S)J 

QGAS (305) = K(305) [ T(605) - T(305) ] 


and 


Q INS (305) = K(105) [ T(525) - T(305) ] 

In a similar fashion the conductances are determined around the model to 
the cooling fin for nodes 304, 303, 302, 301, 201, 202, 203, 204, 401, 

402, 403, 404, 405, 406, 413. Likewise the conductances are determined 
for the other half of the model, nodes 307, 308, 309, 310, 311, 312, 101, 
102, 103, 104, 412, 411, 410, 409, 408, 407, 413. Finally the conductance 
between nodes 413 and 414 is determined. 


6. The program then makes a comparison of the upgraded and input conductances 
by calculating the ratios K , (i)/K(i) where 

K 1 (i) = upgraded conductance determined using experimental data 

K(i) = input conductance determined from kA /Ax 

Figure 11-52 shows the first comparison of the upgraded conductances with the 
input values for the five test runs. This comparison was used to gain the 
insight needed to upgrade the expanded thermal math model. Several observa- 
tions can be made from Figure 11-52: 

1. The conductance ratio is low for the conductors between 

a. The inner cylinder and the end plates, nodes 312 to 101 and 301 to 

201 . 

This discrepancy is explained by the comparison of the input to the 
actual geometry of the end plate to inner cylinder connection. 
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Figure 11-52: INITIAL UPGRADING RUN RESUL TS 
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Figure 11-53 shows this comparison along with an approximate cal- 
culation of the expected conductance ratio. This ratio is in agree- 
ment with that calculated from the experimental data, 

b. The end plates and the outer cylinder, nodes 104 to 412 and 204 to 401. 

This discrepancy is also explained by the difference in the input and 
actual geometry of this connection. Figure 11-54 shows a comparison 
of the input and actual geometry along with an approximate calculation 
of the conductance ratio. This ratio is in fair agreement with that 
calculated from the experimental data. 

c. The cooling fin and the cooling water coils, nodes 413 to 414. 

This discrepancy is apparently due to the input value assuming an 
infinite conductance from the end of the fin to the cooling coils, 

A more reasonable approach , see Figure 11-55 , gives a conductance 
ratio in agreement to that calculated from the experimental data. 

2. The conductance ratio is high for the conductors between 

a. The heater section and adjacent nodes, nodes 306 to 305 and 307 to 
308. 

IcA 

The discrepancy is partially explained by using ( — ) to calculate 
the input conductance for the heated section. The conduction of the 
heater core may also contribute to the discrepancy. 

b . The inner cylinder nodes 311 and 312 

This difference is probably due to the temperature difference around 
the inner cylinder since only one thermocouple is available at node 
312 and the temperature at node 311 is taken as the average of the 
two thermocouple readings at that position. 


413 



INPUT ACTUAL 

GEOMETRY L GEOMETRY 


D180-15048-1 



414 





INPUT GEOMETRY ACTUAL GEOMETRY 


D180-15048-1 





415 


Figure 11-54: END PLA TE TO OUTER CYLINDER CONNECTION 



INPUT actual 

GEOMETRY GFOMFT 


D180- 150 48-1 




416 


Figure 11-55: COOLING FIN GEOMETRY 




D180-15048-1 


3. There is greater scatter in the data on one half of the model than the 
other half. 

The side with the larger scatter has the thermocouple and 
heater leads routed through it and the scatter is probably 
related to the heat leak not being calculated properly. 

4. There is a tendency for the conductance ratio to increase along the outer 
cylinder from the end plate to the fin especially for nodes 412 to 407. 

This may be caused by too much radiation heat transfer in the 
thermal math model. The script F values were calculated 
for an emissivity of 0.88 whereas the measured sample value 
was 0.83 before baking at 400°F and 0.88 after baking. 

5. The conductance ratio is generally less than unity. 

Possibly the input value of thermal conductivity for the 
stainless steel is somewhat larger. 

6. The data scatter appears to be systematic with heat input, especially for 
nodes 201 to 204 and 401 to 405* 

This is probably due to the heat leak not being calculated 
properly. 

7. The scatter in the data is largest in regions away from the heater section 
and cooling fin. 

These are the regions in which the conductive heat trans- 
fer is small and minor error elsewhere are amplified in 
these regions. 
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Table 11-13 outlines the 25 thermal math model upgrading runs. These runs 
investigated the effects of changes in the following aspects of the math model: 

o Radiation exchange factors (emissivity) 

o Air supply and return tube effective conductance to ambient 
o Heater lead routing and conductance 

o Foam insulation conductance 

o Heat leak due to vacuum line penetrations into foam 
o Thermal conductivity of phenolic air distribution plugs 
o Conductance between heater and inner cylinder wall 
o Conductivity of stainless steel 

The second upgrading run reduced the emissivity from 0.88 to about 0.865 
(assuming F ^F 0 = ( The results of this run are shown in Figure 11-56. 

The emissivity change corrected the tendency of the conductance ratio to 
increase along the outer cylinder toward the fin. However, the data scatter 
increased for the end plates and outer cylinder. A more detailed evaluation was 
then made of the heat leaks and the conductances of several heat flow paths were 
changed for run 3. The results of this run, shown in Figure 11-57, show some 
reduction in the data scatter. Runs 4-9 investigate the effects of changes in 
thermal conductivity of stainless steel and phenolic, heater lead routing and 
increased conductance of the air lines. The changes for the phenolic had little 
effect on the results. The increase in thermal conductivity for steel slightly 
increased the scatter. The routing of the heater leads to the outer air 
distribution plug and the increase of the air tube conductances somewhat reduced 
the data scatter. 

Runs 10 and 11 compared the extreme cases of the heater leads having no inter- 
action with the model between the heater and ambient and of the heater leads 
having no direct interaction with ambient. Connecting the heat leads directly 
to ambient, run 10, resulted in reduced data scatter for the top end plate, 
nodes 201-204, and upper half of the inner cylinder, nodes 401-406. However, 
the resulting conductance values for the top end plate appeared unrealistically 

low (k , ,/k. » 0.7). Disconnecting the heater leads from ambient arid 

upgraded input ° 

connecting them directly to node 520, run 11, increased the conductance values, 
but also increased the data scatter. Run 12 was made in an attempt to reduce 
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Script F r s changed by factor of 0.933 from 



leads reconnected to ambient K(284) = 0,010; added conductance 
um lines to insulation, K(285) = 0.010. 
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Changed heater to cylinder conductance, K (106) & K(107) , to 0,20 
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Figure 11-56: UPGRADING RESUL TS WITH REDUCED EMISSIVITY 
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Figure 11-57: UPGRADING RESUL TS WITH REVISED HEA T LEAK PA THS 
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the scatter by increasing the emissivity. This attempt was only partially 
successful and the results inconclusive » 

At this point a complete check was made of the TMM to see if any errors were 
present. It was discovered that the radiation exchange factors (FVs) were not 
quite symmetric. This occurred due to some redundant calculations being made 
for these factors with the results not being quite equal, i.e., the calculated 
values for A were not quite equal to those for A Subsequently these 

radiation exchange factors were changed to be symmetric with the AF value taken 
as the average of the two calculated values. The results of run 13, a repeat 
of run 12 with the revisions to the radiator cards, showed the effects of these 
changes to be very minor. 

The results for most of these remaining runs are shown in Figure 11-58 and 

11-59 which give the average conductance ratio, k , k, . and the range 

& ° upgraded input 1 

of variation for this ratio, the difference between maximum and minimum divided 
by the average value. 

Figure 11-58 shows the effects of changing the penetration conductors K(284) 
(heater leads) and K(285) (vacuum lines) , the conductors between the heater 
and the wall K(106, 107) and one case of emissivity change. The conductance 
of the insulating foam remained at it-s original input value. 

Figure 11-59 shows the effects of changing the emissivity, the foam conductance 
and the conductors between the heater and the wall . 

Inspection of Figures 11-58 and 11-59 indicates that the upgraded conductance 
ratio is quite sensitive to changes in the TMM which are within the range of 
uncertainty. The range of variation in the conductance ratio is not so 
sensitive to these changes, however; the individual test run deviation from the 
average ratio is strongly affected by some of the changes in the TMM and 
different test runs give the extremum variations for a given conductor. 
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Figure 11-58: UPGRADING RESUL TS FOR TEST RUNS 10, 14. 15, 16, AND 17 
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Figure li-59: UPGRADING RESUL TS FOR TEST RUNS TO, 14. 15, 16, AND 17 
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Consideration was given to the detailed results of the upgrading runs in an 
attempt to discover cause and effect relationships between the results. In 
many cases apparent cause and effect relationships for the scatter at the top 
end of the model have the opposite effect at the bottom end of the model. 

The following aspects of the model may be responsible for limiting the degree 
to which the thermal math model could be upgraded: 

o Paint emissivity - The emissivity of the paint may be temperature 
dependent or the paint may not be uniform over the surfaces. 

o Foam Insulation - The foam insulation may not be uniform in properties. 
(The insulation was formed by pouring several batches of foam.) 

o Conductance between heater and cylinder - The conductance between the 

heater and cylinder may be nonuniform and is probably influenced by the 
"joint conductance" between the RTV and cylinder wall. The "joint 
conductance" is strongly dependent on the pressure at the joint and 
consequently strongly dependent on the heater temperature. 


II. 5. 7. 2 Free Convection Test Series Correlations 


II. 5. 7. 2.1 1/4 Scale Model Nusselt Number Correlation 

The 1/4 scale model data for the free convection test series are listed in 
Appendix II-A.1.2. Typical temperature distributions are shown in Figure 11-60 
for different pressures at a given heating rate. 

Application of the Nusselt number preservation scaling technique to the free 
convection tests requires the Nusselt number to be correlated with the Grashof 
number. Appropriate definitions of the Nusselt and Grashof numbers must be 
made in order to achieve a correlation between them. 
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Figure 1 1 -60: TYPICA L FREE CONVECTION TEST TEMPERA TURE DISTRIBUTIONS 
1/4 SCALE MODEL - HIGH HE A TING RA TE 
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The critical factor in these definitions is the temperature differences that 
are used. The Grashof number is generally based on the characteristic tempera- 
ture difference between the wall and the gas and the local Nusselt number is 
generally based on the local temperature difference between the wall and the 
gas. In the model tests there were no measurements of the gas temperature 
distribution, consequently either the gas temperature must be calculated in 
some manner or other temperature differences must be used. It should be noted 
that measurement of the gas temperature distribution would add a large amount of 
instrumentation and would probably result in complicating the problem. 

Two subroutines were developed for the expanded thermal math model to analyze 
the 1/4 Scale Model "Free Convection Test Series" data. One subroutine was 
used to input the upgraded conductance values and the other used to analyze the 
convective heat transfer in the model. The upgraded conductances and other 
pertinent parameters from upgrading run number 20 were used in the data analysis. 
The data analysis subroutine makes a thermal balance at each node and determines 
the convective heat transfer rate. The Nusselt number for each node and the 
Grashof number for the test run are then calculated. Several definitions for the 
Nusselt and Grashof numbers were tried before an adequate correlation was 
achieved . 


The first correlation attempt based the Nusselt number on the difference between 
the node and average gas temperature, i.e.. 


where 



Nu. = Nusselt number at node i 
i 1 

k = gas thermal conductivity at T 

A. = area of node i 

x 

Q = heat transfer to gas from node i 
g 

T = temperature of node i 

T = average gas temperatures 

g 

L = characteristic length (taken as 6 inches for the 1/4 scale model) 
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Using this definition of the Nusselt number and noting that the net heat 
transfer to the gas must be zero, the average gas temperature can be written 
as 



The subroutine used an iterative process to solve these equations for the 
average gas temperature and the local Nusselt numbers. The Grashof number 
was based on the rms temperature difference between the wall and the gas with 
an area weighting factor, i.e., 





The results of this analysis, using the average gas temperature, failed to 
produce a correlation between the Nusselt and Grashof numbers . The calculated 
thermal balances, however, were useful in further correlation attempts. 

The calculated overall thermal balance, using the upgraded conductance values 
from upgrading run No. 20, is shown in Table 11-14 for each of the "Free 
Convection Test Series" runs. An overall thermal balance is achieved to within 
3.5 percent of the heat input for all of the runs. There is, however, a 
calculated net heat transfer to the gas for all of the runs which is about 8 
percent of the heat input for the worse case. This indicates further upgrading 
of the TMM may be required. It is interesting to note that up to 70 percent 
of the heat input is transferred to the gas. 

Since the use of the average gas temperature did not allow a correlation of the 
data the detailed computer results were used to search for a new analysis 
scheme. The proper gas temperature to use in the calculation of convective heat 
, transfer from the common equation Q = hA (T^ - T ) is known at the two locations 
on the model where the heat transfer to the gas is zero. At these points 
T^ = Tg. Having these two gas temperatures it is possible to calculated the gas 
temperature variations around the model from the calculated heat transfer to 
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TABLE 11-14 CALCULATED THERMAL BALANCE FOR 1/4 SCALE MODEL FREE CONVECTION RUNS 
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the gas, assuming a uniform gas circulation, by integrating the heat transfer to 
the gas between the two positions, i.e.. 


" 7 ^ 00 = Tj 60 6(, ) — ~C* (rfpij 




£'<*, <** 

A hand calculation was made for test run 20.1, using the calculated thermal 
balance, to determine the gas temperature variation around the model. Figure 
11-61 shows the model and calculated gas temperature as well as the heat 
transfer to the gas at each affected node. The Nusselt number, based on the 
local temperature difference between the gas and the wall, could probably be 
correlated with the Grashof number. However, it would be difficult to use this 
correlation in the full scale model analysis. Consequently a simpler correla- 
tion technique was sought. 


If the temperature profiles for all of the test runs are similar to those 
shown in Figure 11-61, then the Nusselt number could be based on a characteristic 
temperature difference and a correlation with the Grashof number should be 
obtained. For similar temperature profiles the temperature crossover points 
should remain at the same position. A calculation of these crossover points 
showed that their positions were nearly the same for all of the test runs. 

Table 11-15 lists these positions for each run. The slight variations of 
position may be due to inaccuracies in the TMM. 

A direct correlation of the local heat transfer to the gas as a function of 
Grashof number was tried but was not successful. However, the heat transfer to 
the gas divided by a characteristic temperature difference showed promise of 
correlating with the Grashof number based on this temperature difference. 


Consequently the data analysis scheme was revised to base both the Nusselt and 
Grashof numbers on a characteristic temperature difference, i.e.. 



and 


<V4 


Gr = 


ptgeLf/xr, 

A* 
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Figure 11-61: DA TA ANAL YSES RESUL TS FOR TEST RUN 20. 1 
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TABLE 11-15 

CALCULATED POSITIONS WHERE HEAT TRANSFER TO GAS IS ZERO 


Outer Cylinder Wall - 
Run Number Inches From Top 


Inner Cylinder Wall - 
Inches From Top 


18 

4.15 

1.48 

19 

4.20 

1.45 

21 

4.15 

1.50 

20.1 

4.35 

1.46 

37.1 

4.24 

1.42 

38.1 

4,43 

1.37 

39 

4.42 

1.31 

22 

4 .22 

1.51 

23 

4.34 

1.48 

34 

4.40 

1.38 

35.1 

4.46 

1.40 

36 

4.45 

1.32 

25 

4.49 

1.51 

24 

4.46 

1.49 

31.2 

4.37 

1.46 

32.1 

4.48 

1.42 

33 

4.44 

1.34 

26.1 

4.51 

1.57 

27.1 

4.44 

1.54 

28 

4.42 

1.50 

29.1 

4.48 

1.44 

30 

4.51 

1.40 
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The characteristic temperature difference was defined as one half of the 
difference between the average inner and outer cylinder temperatures, 

AT„ = i - 72au.) 

The results of this analysis produced a correlation between the Nusselt and 
Grashof numbers. These results are listed in Appendix II-B . Since the Nusselt 
number is based on a characteristic temperature difference negative values are 
calculated for nodes which are heated by the gas and positive values for those 
cooled by the gas. The results were correlated by plotting the Nusselt number 
versus the Grashof number for each inner cylinder node, pairs of outer cylinder 
nodes and the four nodes (area averaged) of each end plate. As expected the data 
scatter for the nodes having large convective heat transfer was much less than 
that for those having a small convective heat transfer. Figure 11-62 shows 
typical results for nodes having relatively large convection effects. The 
curves for the nodes are quite linear, however, their slopes differ. Figure 
11-63 shows typical results for nodes having small convection effects. The 
data scatter is relatively large and the correlation is marginal. 


Since the variation of the heating rate at a given pressure does not produce 
large variations in the Grashof number the data were averaged at each pressure 
and the results plotted for all the nodes. Figure 11-64 shows the correlation 
for the inner cylinder nodes and Figure 11-65 shows the correlation for the outer 
cylinder and end plate nodes. 

The Nusselt number variation around the model is shown for three Grashof numbers 
in Figure 11-66. From this Figure it is apparent that the regions where the 
data scatter is large have negligible convection effects. Consequently the 
data correlation should be adequate for predicting the full scale model 
performance. 

The Nusselt number preservation scaling technique investigation was terminated 
at this point. There was not adequate time to pursue this investigation to its 
completion. However, the work completed does allow an assessment of this 
scaling technique . 
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Figure 11-63: TYPICAL NUSSEL T NUMBER CORRELA TIONS FOR REGIONS 
WITH SMALL CONVECTION EFFECTS 
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Figure 11-66: NUSSELT NUMBER VAR! A TION AROUND MODEL 
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The Nusselt number correlation developed for the free convection test series 
using the 1/4 scale model can be used with an expanded thermal math model of 
the full scale model to predict its thermal performance. This expanded math 
model would make use of the information gained from upgrading the thermal math 
model for the 1/4 scale model. In particular the conductances at the weld joints 
would be reevaluated. 

It should be possible to develop a Nusselt number correlation for the forced 
convection test series using the 1/4 scale model data. This correlation is 
expected to be more complex and more difficult to develop than that for free 
convection. The amount of time and effort required to develop this correlation 
cannot be accurately predicted. 

II. 5. 7. 2. 2 Correlation Between Full Scale and 1/4 Scale Models 

The results of the 1/4 scale model free convection tests made at a pressure of 
1/2 atmosphere were compared directly with the full scale model free convection 
tests (the full scale model free convection test series data are listed in 
Appendix II-A.2.1). This comparision is shown in Figure 11-67. The full scale 
model temperatures are somewhat greater than those for the 1/4 scale model. This 
temperature difference is greater at the higher heating rates. The temperature 
difference between the full scale and 1/4 scale models is shown for the highest 
and lowest heating rate cases in Figure 11-68. 

The temperature measurements for the full scale model heater section (nodes 
306 and 307) appear to be somewhat erroneous. This was expected because of 
the earlier problems with these thermocouple installations. Apart from these 
nodes the maximum temperature difference for the highest heating rate is about 
13°F and that for the lowest heating rate 3.5°F. These temperature differences 
are probably caused by slightly different emissivities in 1/4 scale and full 
scale model. The paint sample measurements taken after "baking 11 gave an 
emissivity of 0.83 for the aluminum base and 0.88 for the stainless steel base. 

It is estimated that this difference could result in a higher full scale model 
temperature by as much as 20°F at the heater section for the highest heating 
rate case. Other differences between the 1/4 scale and full scale models do not 
explain the measured temperature differences. The differences in conduction 
between the models tends to lower the full scale model temperatures relative 
to those of the 1/4 scale model. This is the reverse of what is observed. 


440 



D180-15048-1 



441 


Figure / 1-67: THERMAL SCALE MODELING CORRELA TION BETWEEN 1/4 SCALE 
AND FULL SCALE MODELS - FREE CONVECTION TESTS 
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Figure 11-68: TEMPERA TURE DIFFERENCE BETWEEN FULL SCA LE AND 
1/4 SCALE MODELS - FREE CONVECTION TESTS 
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The excess heat leak (due to the heater leads) of the 1/4 scale model is not 
large enough to account for the temperature differences. The maximum effect 
of this heat leak is estimated to be 2-3°F. The temperature differences 
between the models does not appear to be caused by differences in the free 
convection heat transfer. Based on the Nusselt number versus Grashof number 
correlation, the average convective heat transfer in the full scale model 
should be about 4 percent greater than that in the 1/4 scale model. 

The correlation between the models shows that the radiation and conduction as 
well as the free convection heat transfer processes can be preserved in a 
thermal scale model. 

II. 5. 7. 3 Forced Convection Test Series Correlations 

The 1/4 scale model data for the forced convection tests are listed in Appendix 
II-A.2.* The full scale model tests were made at atmospheric pressure for 
nominal Reynolds numbers of 300, 1200 and 4800 and nominal heating rates of 
187, 224 and 294 watts. The high heating rate case was repeated for the three 
flow rates with the flow direction reversed. The 1/4 scale model was tested 
at various pressures and heating rates for Reynolds numbers equal to those of the 
full scale model. These tests were to be used to develop a Nusselt number 
correlation for the Nusselt number preservation sealing technique. However, 
as discussed previously no correlation was attempted since adequate time to 
accomplish this task was unavailable. The 1/4 scale model was also tested under 
mass flux and heat transfer coefficient preservation conditions. These tests 
were used for direct correlations between the 1/4 scale and full scale models. 

The correlations between the 1/4 scale and full scale models, at the high 
heating rate, are shown in Figures 11-69, 11-70 and 11-71 for the high, 
intermediate and low flow rate cases respectively. The 1/4 scale model data 

•* The gas outlet temperature measurements for the 1/4 scale model are 
inaccurate especially at the low flow rates. 
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are shown for both the mass flux and the heat transfer coefficient preservation 
techniques. Figures 11-72, 11-73 and 11-74 show the same correlations with the 
flow direction through the model reversed (the flow is from top to bottom 
instead of bottom to top) . These correlations show that the preservation of 
mass flux results in 1/4 scale model temperatures generally lower than the full 
scale model while preservation of the heat transfer coefficient results in 
higher model temperatures. These results are in qualitative agreement with the 
approximate analysis presented in Section II. 4. 3, Figure 11-75 shows the 
results of this analysis using parameters applicable to the model configuration. 
Since the inlet gas temperature is less than the reference temperature (tempera- 
ture without gas flow) these results show that preservation of mass flux results 
in lower 1/4 scale model temperatures and heat transfer coefficient preservation 
results in higher temperatures . 

The temperature distribution correlation between the 1/4 scale and full scale 
models, shown in Figures 11-69 through 11-74, can be qualitatively explained 
as follows: 

At the high flow rate forced convection dominates and the temperature 
correlation for the heat transfer coefficient preservation is directly 
related to the air temperature change through the model. The largest 
temperatures differences occur in the air outlet region and these 
differences closely correspond to the air outlet temperature differences 
between the 1/4 scale and full scale model. At the low flow rate free 
convection dominates and the temperature differences are related to 
the differences in the air cooling rate. The average air temperature 
in the models is fixed primarily by the free convection effects and 
should be nearly equal in 1/4 scale and full scale models. Consequently 
the relatively smaller mass flow rate in the 1/4 scale model results 
in less air cooling and more heat conduction to the cooling fin. This 
results in generally higher temperatures in the 1/4 scale model. At 
the intermediate flow rate combined convection is important and the 
correlation is not as good as that at the high and low flow rates. 
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Preservation of mass flux results in good correlation for the end 
plates and outer cylinder temperatures. Since the full scale model 
end plates and outer cylinder temperatures are not ,far different from 
the air temperature at the higher flow rates, the increased convective 
heat transfer in these regions tends to bring these temperatures even 
closer to the air temperature. This tendency is offset somewhat by 
the increased air cooling of the inner cylinder which results in less 
heat transfer from the inner to outer cylinder. The net effect then 
is for the mass flux preservation technique to result in good correla- 
tion for the end plates and outer cylinder and poorer correlation for 
the inner cylinder which excessivley cooled by the air. As expected 
the correlation is best for the low flow rate case where free convec- 
tion dominates. The correlation for the inner cylinder improves as 
the flow rate is decreased. 

The correlations between the 1/4 scale and full scale models for the intermediate 
heating rate and presented in Figure 11-76 which shows the temperature differ- 
ences between the 1/4 scale and full scale model at each of the flow rates 
for both scaling techniques. This Figure shows that the two techniques do not 
quite bound the full scale model temperatures. The higher convective heat 
transfer coefficient for mass flux preservation keeps the outer cylinder warmer 
in the regions near the cooling fin. 

The forced convection test series data were compared with the approximate 
analysis presented in Section 11.4.3 by averaging the end plate, outer cylinder 
and inner cylinder node temperature. This comparison was made for the three 
flow rates and the three heating rates. The reference temperature for a partic- 
ular heating rate was taken as the average temperature of the full scale model 
under free convection conditions at the same heating rate. Figure 11-77 shows 
the nondimens ional temperature difference between model and prototype as a func- 
tion of Reynolds number for the three heating rate cases. Also shown in Figure 

g 

11-77 are the results, calculated for a prototype Grashof number of 10 , for a 
1/5 scale spacecraft (see Section II. 4. 3) as well as those for the 1/4 scale 
model. Even though the analysis is a poor representation of the model the 
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calculated results are in fair agreement with the data, however, the maximum 
temperature differences for heat transfer coefficient preservation do not coin- 
cide. Since the approximate analysis is more representative of a manned space- 
craft configuration than of the model configuration the calculated results should 
give a good indication of the temperature differences to be expected in thermal 

scale modeling a manned spacecraft* The results shown in Figure 11-77 indicate 

4 

that, at a typical manned spacecraft Reynolds number of 10 , better thermal 
similitude would be achieved in a 1/5 scale spacecraft than that achieved in the 
1/4 scale model. 
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II. 6 CONCLUSIONS 


Adequate thermal similitude may be achieved in thermal scale modeling of 
radiation-conduction-convection systems by using the compromised scaling 
techniques of either mass flux or heat transfer coefficient preservation. The 
best thermal similitude is achieved with mass flux preservation when free con- 
vection effects dominate and with heat transfer coefficient preservation when 
forced convection effects dominate. The degree of thermal similitude achieved 
with either technique depends on the system being modeled. This is illustrated 
by the approximate analysis presented in Section II. 4. 3 and was verified in the 
experimental investigation* 

The heat transfer coefficient preservation scaling technique should give better 
thermal similitude than mass flux preservation in thermal scale modeling applica- 
tions for manned spacecraft. The approximate analysis of Section II. 4. 3 indicates 
that very good thermal similitude (both transient and steady state) can be 
achieved with a 1/5 scale model spacecraft using heat transfer coefficient 

preservation. This analysis indicates that, for a typical manned spacecraft 

4 

Reynolds number of 10 , the temperature agreement between 1/5 scale model and 
prototype spacecraft will be within about 5°F and that the transient response 
time will be scaled to within about 5 percent. 

The major problem involved in using heat transfer coefficient preservation is 
determining the proper relationship between the Nusselt and Reynolds number. 
Transitions between laminar and turbulent flow may also cause problems. However, 
In the experimental investigation flat plate laminar flow theory was used for the 
relationship between Nusselt and Reynolds numbers with good results even through 
the inlet flow was highly turbulent and the flow through the model was in the 
laminar flow regime. The experimental investigation represents a more severe 
test of thermal scale modeling involving convection than does thermal scale 
modeling of a manned spacecraft. These more severe test conditions were large 
conduction and radiation effects and the high inlet gas velocities. Consequently 
the use of heat transfer coefficient preservation should not present any major 
problems in thermal scale modeling of manned spacecraft. 
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The Nusselt number preservation scaling technique is a workable technique that, 
if carried out properly, will allow prototype performance to be predicted using 
the Nusselt number correlations developed from the thermal scale model tests. 

Much of the detailed work reported in this document indicates the methods required 
to implement this scaling technique. Thermal scale modeling using this technique 
must be approached as an experimental research program rather than as an engineering 
test program. As is the case with most research investigations involving convective 
heat transfer great care is required throughout the investigation to achieve an 
accurate determination of the heat transfer coefficients. Once the heat transfer 
coefficients are accurately determined, the use of this techinque assumes that a 
valid correlation for the Nusselt number can be found. This correlation of the 
Nusselt number with the system parameters requires the investigator to be thor- 
oughly versed in convective heat transfer. The experimental investigation carried 
the Nusselt number preservation scaling technique to the point of achieving a 
correlation for the 1/4 scale model free convection data. There was insufficient 
time to carry this investigation to its completion. It is estimated that at least 
an order of magnitude greater effort is required to use the Nusselt number preser- 
vation technique than the other techniques. For some critical problem areas it may 
be necessary to use Nusselt number preservation to accurately determine the 
convection effects. However, before such a task is undertaken the requirements 
of using this technique should be well understood since any attempts to do a 
"quick and dirty" job will result in a wasted effort. 
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II-A.l 1/4 SCALE MODEL DATA 

II-A.1.1 RADIATION-CONDUCTION TEST SERIES 
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RUN NUMBER 39.0 VOLTAGE 33.67 VOLTS CURRENT 0*2597 AMPS POwER 8.747 V* A T T S PRESSURE 7. 930 ATM FLOW 0.0 LB/ MIN 
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RUN NUMBER 50.0 VOLTAGE 48.87 VOLTS CURRENT 0.3760 AMPS POWER 18.376 WATTS PRESSURE 1.078 ATM FLOW 0.28500 LB/M 1 N 
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RUN NUMBER 62.0 VOLTAGE 48. 93 VOLTS CURRENT 0.3757 AMPS POWER 13.385 WATTS PRESSURE 0.097 ATM FLOW 0.U74S0 LB /MI N 
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III.l INTRODUCTION 

The spacecraft docking thermal interface may be critical when a controlled 
temperature spacecraft is being docked to a spacecraft without temperature 
control. For example the heat transfer from the controlled temperature 
spacecraft to a cold spacecraft during docking may result in local moisture 
condensation on the inside wall of the controlled spacecraft. Disregarding 
the possibility of convective heat transfer between the docked spacecraft, 
the spacecraft docking thermal interface involves both radiative and conduc- 
tive heat transfer. The thermal interface may be separated into a "local 11 
interface and a "distant" interface. The "local" interface involves both 
radiation and conduction, whereas; the "distant" interface involves only 
radiation heat transfer. 

Thermal scale modeling application to the "distant" interface is within the 
present state-of-the-art. This modeling application can be used to simulate 
the radiant interchange and blockage effects for complex multiple docking 
configurations. Thermal scale modeling application to the "local" thermal 
interface may be more difficult due to the possible complexities in the dock- 
ing mechanisms . In particular the contact conductance at the docking inter- 
face may be difficult to simulate in a thermal scale model. 

This report presents the results of an investigation of transient thermal 
scale modeling applied to the spacecraft docking thermal interface. The 
thermal scale modeling criteria and scaling techniques applicable to the 
docking interface are described. An investigation of a typical AAP docking 
configuration was made and the details of the spacecraft docking thermal 
interface were defined. This definition includes the conduction and radia- 
tion heat transfer processes with particular emphasis on the contact con- 
ductance at the mating flanges. 
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Ill % 2 THERMAL SCALE MODELING CRITERIA 


The thermal scale modeling criteria for the spacecraft docking thermal 
interface are the same as those for any radiation-conduction system with 
the additional complexity of the contact conductance scaling criteria. The 
radiation-conduction scaling criteria are given in Part II of this report. 
The contact conductance scaling criteria are easily derived by relating 
the contact conductance to the heating rate per unit area, 

h AT = q (1) 

where 

h = contact conductance coefficient 
c 

AT = temperature difference across joint 
q = heating rate per unit area 

The scaling criteria for q is that (q/aT^) remain invarient , where T q 
is the characteristic temperature. Consequently the scaling criteria for 
contact conductance is 



where the superscript * refers to the model to prototype parameter ratio. 

For the temperature preservation scaling technique 

h * = 1 (3) 

c 

and for material preservation 

h * - (T V - (lV 1 <4> 

c o 

where L is the characteristic length. 
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Rice (reference 1) has shown that the contact conductance can be represented 
by an equation of the form 



(5) 


where 

c = constant 

k = mean thermal conductivity 
m = mean rms slope of surface roughness 
a = mean rms height of surface roughness 
P = pressure at interface 
H = micorhardness of material 
a = constant 


The microhardness H may be approximated, in terms of yield strength S, by 
H ^ S/3 (6) 

and by plotting experimental data on log-log coordinates Rice obtained 

4- O.M S(§)(%)"~ 

This equation holds for ideal surfaces in contact, i.e., surfaces which 
do not have waves or distortions. Surfaces of T 0 f ring groove quality are 
examples of the type surface for which the contact conductance is given by 
equation (7). 

The conditions required for meeting the contact conductance scaling criteria 
may be found from equations (3), (4) and (7). 

The criteria for temperature preservation scaling then becomes 
or in terms of pressure 



512 



D180-15048-1 


If the mating surfaces of both the prototype and model are machined to the 
same roughness then the scaling requirement for temperature preservation 
becomes 


* 

P 


* 

S 


(kV 1 * 18 


( 10 ) 


The required pressure at the mating surfaces for temperature preservation 
scaling is thus a function of the yield strength and thermal conductivity. 


The criteria for material preservation scaling can be written using equations 
(4) and (7) as 

(&) *(p«) 0 ' 9 B =(l*)'' (11 > 

or in terms of pressure at the interface, assuming identical surface 
conditions for model and prototype. 


* 

P ~ 


* 


(L ) 


-1.18 


( 12 ) 


The material preservation scaling technique thus requires the interface 
pressure to be increased to model the contact conductance. 
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III. 3 DOCKING INTERFACE DEFINITION 


The docking interface definition was based on the Apollo Command Module (CM) 
docking system* In this study the thermal effects of the docking probe and 
drogue were considered insignificant and emphasis placed on the heat transfer 
through the basic docking ring structures. In the docking of the Apollo CM 
to another spacecraft the Apollo docking ring is held to the mating spacer- 
craft docking ring by twelve spring loaded latches. The docking ring is a 
35 inch diameter tube with flanges at both ends. There are two 0-ring grooves 
in the mating flange and one O-ring groove in the flange which is held per- 
manently to the CM with metal to metal contact along its surface. 

The heat transfer between the docked spacecraft occurs primarily by conduc- 
tion through the docking ring structures and radiation through the docking 
tunnel. The most complex heat transfer process occurs at the mating flanges 
where contact conductance effects may be significant. These contact con- 
ductance effects depend on the details of the docking latches, in particular 
on the amount of docking force generated by the latching system. In the 
accomplishment of this task the details of the docking latches were unavail- 
able until late in the program. Consequently the docking interface was 
defined, except for the latches, and thermal analyses conducted for various 
values of contact conductance at the mating surfaces. When the details of 
the docking latches became available their effect on the interface was also 
determined. 

Ill .3.1 Thermal Analyses 

A thermal math model was developed for the docking interface between the CM 
and the Subsystem Test Bed (STB) . Since the conduction through the Apollo 
docking ring is a major heat transfer mode for the thermal interface, a 
detailed analysis was performed to determine the total thermal conductance 
of the CM docking ring. The results of this analysis was then used to 
design a thermally equivalent docking ring that replaces the complex CM 
docking ring cross section with a simple cross section. This thermally 
equivalent CM docking ring cross section is shown in Figure III-l . The 
overall dimensions are the same as the actual docking ring. 
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FIG °E 111-1 THERMALLY EQUIVALENT CM DOCKING RING 


515 


D180-15048-1 


The thermal math model nodes used for the CM/ STB thermal interface is shown 
in Figure III-2 . This math model includes convection from the STB atmosphere 
to the STB pressure shell and docking hatch; conduction along the STB pres- 
sure shell to the docking tunnel; conduction along the tunnel wall; variable 
contact conductance across the mating interface and conduction along the CM 
docking ring to the CM structure. The radiation heat transfer inside the 
docking tunnel was also included in the math model. The radiative inter- 
change factors (script - F values) for the tunnel, nodes were determined from 
The Boeing Company’s Thermal Radiative Interchange Factor Computer Program 
(reference 2) . The STB tunnel surfaces were assumed to have an emissivity of 
0,2 and the CM thermal hatch surface an emissivity of 0.4, The radiation 
exchange between the mating flanges was also included for cases in which there 
is no metal to metal contact. 

The conductive heat transfer across the mating flanges was considered for two 
cases. In one case the 0-rings are not completely compressed and there is no 
metal to metal contact at the interface. In this case the conduction across 
the interface depends on the degree to which the 0-rings are compressed. As 
the distance between the mating flanges is decreased the contact area between 
the surfaces and the 0-rings is increased due to the flattening of the 0-rings 
This increase in contact area and reduction in conduction path length was 
included in the determination of the 0-ring conductance versus the CM/ STB 
interface separation shown in Figure III-3 . The band of values shown is due 
to the uncertainties in exact sizes and materials used in the docking system. 

The other case considered occurs when the 0-rings are completely compressed 
and the mating flanges contact each other. In this case there is metal to 
metal contact and the dominant mode of heat transfer is by contact conductance 
This contact conductance is very sensitive to the contact pressure exerted at 
the interface (see equation 7) . 

Thermal analyses of the CM/STB docking interface were made using the Boeing 
Engineering Thermal Analyzer (BETA) computer program (reference 3) . These 
analyses assumed the gas atmosphere in the STB to be maintained at 530°R and 
the CM structure to be at 490°R. Three convective heat transfer coefficient 


516 




D180- 15048-1 


FIGURE 111-3 O-RING CONDUCTANCE 
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values (0*05> 0.5 and 1,0 BTU/FT^-HR-°F) were used for the STB atmosphere/ 

wall thermal interface and various values of conductance, ranging from 0 to 
4 

10 BTU/HR-°F, were used for the mating flange interface. Typical steady 

state temperature distributions calculated for the docking tunnel are shown 

2 

in Figure III-4 for a convective heat transfer coefficient of 0.5 BTU/FT - 
HR-°F. This figure shows that there is a substantial temperature drop at 
the mating interface for all cases except for intimate metal to metal con- 
tact where the interface conductance is greater than 100 BTU/HR-°F. The 
heat transfer rate between the STB and CM is shown in Figure III-5 as a 
function of the conductance at the mating interface for the three heat trans- 
fer coef f icients . This figure shows the heat transfer rate to be independent 
of conductance for both high and Low conductance values. This was expected 
since at high values of interface conductance the heat transfer rate is fixed 
by the conduction in the tunnel wall and at low values of interface conduc- 
tance it is fixed by the radiation heat transfer. In the intermediate range 
of contact conductance (from about 1. to 100 BTU/HR-°F) the magnitude of this 
conductance has a significant effect on the spacecraft docking thermal inter- 
face. If the contact conductance is either large or small then its magnitude 
does not have a significant effect on the docking interface. 

III. 3. 2 Conductance At The Mating Interface 

When the details of the docking latches became available it was determined 
that the total conductance for the twelve latches is about 0.5 BTU/HR-°F if 
they are made of stainless steel and about 4.0 BTU/HR-°F if made of aluminum. 
This conductance is important only if there is no metal to metal contact at 
the mating flanges. 

The contact conductance at the mating flanges may be found from equation (7) 
using the following values for the CM/STB docking interface 

k = 79 BTU/FT-HR-°F 
S = 45,000 PSI 

_r 

a - 32 x 10 inches 
m - 1 
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STB GAS TEMPERATURE = 530 DEG. R 
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The contact conductance coefficient at the mating interface is then given as 
h c = 710 p 0 ' 85 BTU/FT 2 -HR-°F (13) 

where P is the interface contact pressure in PS I. 

The contact conductance (h A) at the interface is thus related to the net 

c 

docking force by 

K = 4.9F 0,85 A 0 ' 15 BTU/HR-°F (14) 


where 


K = contact conductance at mating flanges 
F = net docking force lb 

2 

A = contact area of mating surfaces in 

The net docking force is the force exerted by the docking latches minus that 
required to compress the 0-rings (about 1000 lbs) and that required to over- 
come the pressurization of the tunnel (about 5500 lbs) * The total force 
exerted by the docking latches is reported to be about 38,000 lbs. Con- 
sequently the net docking force is about 31,500 lbs. Using this net force 

2 

along with the flange area of 94.5 in in equation (14) gives a contact con- 
ductance of 

K = 6 x 10 4 A a °‘ 15 BTU/HR-°F (15) 

where A^ - fraction of total area in contact. 

Equation (15) shows that even with a small fraction of the flange area in 
metal to metal contact the value of the contact conductance is large enough 
to assume intimate contact. Consequently contact conductance has an 
insignificant effect on the docking interface. 
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III -4 CONCLUSIONS 


Thermal scale modeling of systems involving contact conductance is a problem 
area where further investigation is needed. However the contact conductance 
at the mating flanges is too large to have a significant effect on the Space- 
craft Docking Thermal Interface. Consequently the application of thermal 
scale modeling to this interface presents no special problems and thermal 
similitude may be achieved using existing scale modeling techniques. 
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PART IV 

PRELIMINARY STUDY FOR THERMAL SCALE 
MODELING OF LIQUID LOOP SPACE RADIATORS 
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c 


D 

e 

g 

f 


f 


c.p 


f (Re) 


g 

h 


Gr 


Gz 


k 

L 

l f 

Nu 


Nu(Gz) 

Nu 

q 


p 

p* 


Nomenclature 
specific heat of fluid 
tube diameter 

unit vector in direction of gravity 
friction factor 

constant property friction factor 

friction factor as function of Reynolds number 

acceleration of gravity 

heat transfer coefficient 

Grashof number 

PgS (T - T ) D 3 
w 

2 

y 

Graetz number 

thermal conductivity of fluid 

characteristics length 

radiator fin width (or tube spacing) 

Nusselt number 

hD 

k 

Nusselt number for constant fluid properties 
Nusselt number as function of Graetz number 
Nusselt number for constant wall heat flux 
pressure 

nondimensional pressure 
T 2 

P ... L - £ 

2 

y 
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Pr 


Q 

Re 


t 

T 

T. 

x 

T 


o 



V 

V 


* 


w 

X 


Prandtl number 

PC 

k 

heat rejection rate 

Reynolds number 

pvD _ 4w 

p ttDp 

time 

fluid temperature 

inlet fluid temperature 

characteristic temperature 

radiator fin root (tube wall) temperature 

radiator sink temperature 

fluid velocity 

nondimensional velocity 

v/v 

fluid mass flow rate 

axial distance from tube inlet 


Greek Symbols 

g - coefficient of thermal expansion 

e = emissivity 

= radiator fin effectiveness 

0 ~ nondimensional temperature 

T - T 

= w 

T - T 

w 
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e T 

M 

P 

a 

T 


nondimensional temperature 

T_ 

T 

o 

fluid viscosity 

fluid density 

S tef am-Bol t zmann constant 

nondimensional time 

Kt 

T 2 
P cL 


Subscripts 

m refers to model 

p refers to prototype 

R refers to fin root 

w refers to wall 

Superscript * Refers to model to prototype ratio 

bars denote average values (or evaluation at average fluid temperature) 


Operators 

V 


gradient operator 


e\ a 

1 ax. 

1 


nondimensional gradient operator 
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IV. 1 INTRODUCTION 

Liquid loop space radiators are critical elements in the thermal control of 
manned spacecraft. These radiator systems are quite complex and costly design 
problems have been associated with their development. Thermal scale modeling 
may represent a means of minimizing the development costs involved in the 
design, fabrication and testing of large, complex liquid loop space radiator 
systems . 

A literature search conducted at the start of this study failed to reveal 
any thermal scale modeling studies of liquid loop radiators reported in 
the literature. However, near the end of this study, two spacecraft radiator 
scale modeling investigations were found. In a recent study Colvin and 
Maples, reference 9, have initiated experiments with thermal modeling of 
fluid flow in tubes with a radiation boundary condition. They derived 
"scaling criteria" based on the thermal conduction in the tube wall and 
the convection and radiation boundary conditions for the tube wall. Based 
on these "scaling criteria" they designed and tested prototype and model 
radiator tubes. Their test results for identical Reynolds numbers in model 
and prototype showed a lack of thermal similitude. This lack of similitude 
is related to their neglection of the fluid thermal energy balance in deriving 
"scaling criteria." Fluid energy balance considerations are of primary 
importance in deriving proper scaling criteria. More recently Dietz and 
Fleming, reference 10, have derived scaling criteria and performed an 
analytical investigation of scale modeling a space shuttle radiator panel at 
a 0.4 scale ratio using two scaling techniques. One scaling technique uses 
a fluid change model and the other uses a fluid preservation model. The 
results of their analytical investigation show good agreement between the 
model and prototype temperatures during transient and steady state operation 
for both scale modeling techniques. The scaling techniques described by Dietz 
and Fleming are closely related to those described in this report. 

This report presents the results of a preliminary investigation of thermal 
scale modeling of liquid loop radiators. The characteristics of fluid flow 
in tubes are reviewed and their effects on liquid loop radiator performance 
and tests are discussed. The thermal scale modeling criteria are developed 
and possible scaling techniques discussed. Based on this preliminary study, 
recommendations for further work are made. 
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IV. 2.0 FLUID FLOW IN TUBES (A BRIEF REVIEW) 

Investigations of heat transfer for fluid flow in tubes date back to the 
1870 f s and those of pressure drop date back to even earlier times. These 
investigations are still continuing and accurate predictions for heat 
transfer and pressure drop for all flow conditions are not yet possible 
even though a great amount of literature has been published on the subject. 

Aside from the transitional flow regime, which most designers avoid, the 
major uncertainties in predicting the heat transfer and pressure drop are 
caused by the variable fluid property effects. This is especially true 
for laminar flow where theoretical solutions have been verified for constant 
properties fluids. These variable property effects are less important for 
turbulent flow, however, the theory for turbulent flow is less developed than 
for laminar flow and the emperical correlations are subject to considerable 
uncertainty. 

The variable fluid properties having the largest effect on heat transfer and 
pressure drop are the viscosity and, when in the presence of a gravity field, 
the density. In fluids such as water that have a well defined freezing 
point the variable viscosity effects on the heat transfer coefficient and the 
local pressure gradient are generally small. However, in fluids such as 
oils, which are characterized by a pour point instead of a freezing point, 
the variable viscosity effects can be quite significant. The density 
variations give rise to free convection effects that can be significant for 
all fluids in a gravity field environment. 

IV. 2.1 Governing Equations 

The heat transfer and pressure drop are governed by the energy and momentum 
equations. The energy equation can be written (for constant thermal 
conductivity, negligible viscous dissipation and no internal heat sources) 
as 

pc + pc tr- irr= ^ (i) 
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The momentum equation, Navier-Stokes equation for imcompressible fluid, 
may be written, considering buoyancy as the only body force, as 

?( ~ +&•?&) ~-Vp ~ vxfrvxtr) ( 2 ) 

These equations may be written in nondimens ional form as: 
energy 

dB _ _ .2 - (3) 

■+ ^ ^ ^ ^ & 

momentum 

% 0 + (4) 

These equations are used directly for laminar flow problems, however, for 
turbulent flow, time averaged terms are used and the flow fluctuations are 
accounted for by the addition of turbulent stress terms (Reynolds Stresses) . 

IV .2.2 Laminar Flow 

Numerical solutions have been developed for the energy and momentum equations 
for laminar flow in tubes under various boundary conditions. These solutions 
generally relate to constant property fluids under steady state conditions. 
Tribus and Klein, reference 1, reviewed the solutions for prescribed wall 
temperature boundary conditions. Siegel, Sparrow and Hallman, reference 2, 
present the solution for the uniform wall heat flux boundary condition. The 
radiant heat flux boundary condition case was solved by Chen, reference 3. 

He also relates his numerical solution for the local Nusselt number to 
that for the uniform wall heat flux case to within +2 percent by a simple 
relationship. These results are of particular interest in the design of 
liquid loop space radiators. 

The verification of these numerical solutions with experimental data has been 
limited by the variable fluid property effects on heat transfer and pressure 
drop. These effects have also obscured many experimental data correlation 
attempts. The variable viscosity effects are generally taken into account 
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by using the emperical correlations developed by Sieder and Tate (reference 4). 
Their correlations for laminar flow are: 


heat transfer 




;; \ 0./4 


pressure drop 


f (£) 


- 0.25 


They also correlated turbulent flow data with the same viscosity ratio exponent 
for the Nusselt number correction and with an exponent of -0.14 for the friction 
factor correction. More recently Shannon and Depew, reference 5, developed 
a numerical solution for the variable vicosity effects for laminar flow in 
a tube with uniform wall heat flux. Their solution, which was verified with 
experimental data, shows that the viscosity ratio exponent is a function of 
the Graetz number and is also a weak function of the viscosity ratio itself. 

For fully developed flow the Nusselt number correction is in close agreement 
with the Sieder-Tate correlation for viscosity ratios greater than about 0.1. 
However, the Sieder-Tate correlation for the friction factor does not fully 
account for the variable viscosity effects. The numerical analysis and 
experimental data show that the viscosity ratio exponent is about -0.5 for 
fully developed flow as compared to the Sieder-Tate value of -0.25. 


The free convection effects on heat transfer and pressure drop depend on 
the tube orientation. Vertical tubes lend themselves to theoretical analyses 
because of the symmetry, whereas, the asymmetric effects in horizontal tubes 
make the problem much more complex. Free convection effects in vertical 
tubes have been investigated by Hallman, reference 6, for the uniform wall 
flux case and by Rosen and Hanratty , reference 7, for the uniform wall 
temperature case. These investigations related the free convection effects 
to the magnitude of the Grashof to Reynolds number ratio, Gr/Re . 

Shannon and Depew, references 5 and 8, experimentally investigated the free 
convection effects in horizontal tubes with a uniform heat flux boundary 
condition. The results of these investigations show that the free convection 
effects are independent of the Reynolds number and depend on the magnitude 
of the Rayleigh number (Ra « GrPr) for fully developed flow. Free convection 
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3 

effects become significant for a Rayleigh number between 2 x 10 and 5 
3 

x 10 , The free convection effects increase the Nusselt number by as much 

A 

as a factor of 2 at a Rayleigh number of 2 x 10 . 


In many practical applications both the variable viscosity and free convection 
effects are significant. These combined effects make theoretical analyses 
very difficult and obscure the emperical data correlations. 


IV. 2. 3 Turbulent Flow 


Since the temperature difference between the fluid and tube wall is generally 
small for turbulent flow, the free convection and variable viscosity effects 
are generally small. However the theory for turbulent flow heat transfer 
is less developed that that for laminar flow and there remains some uncertainty 
in the emperical data correlations. The turbulent flow Nusselt number 
correlations generally refer to an average Nusselt number for the entire tube 
length and are expressed in the form 


— n rfl fu / , \ 

NJu = C Re Pr £- ) a(U\ 

i/W ( 710 / 


Various correlations use different values for the parameters in this equation. 
Typical values for the constant C range from .021 to .025; the exponent n on 
the Reynolds number from 0.7 to 0.9; the exponent m on the Prandtl number 
from 0.3 to 0.6; the viscosity ratio exponent 1 for liquids from 0.1 to 0.2; 
and the function of length to diameter ratio g(L/D) is given in various forms. 
Perhaps the most widely used correlation is the Sieder-Tate equation 



Q'OZS Re. 


o.e 



6. 1 4 


( 6 ) 
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IV. 3 LIQUID LOOP SPACE RADIATOR PERFORMANCE 

Consider the simple case of a circular tube attached to a radiator fin with 
one side radiating to a constant temperature sink. The thermal balance 
between the fluid heat losses and the convenction to the radiator fin is 
given by 



and the balance between the convective and radiative heat transfer is 
given by 

rt-kHu (t-t k )= (r^-% 4 ) m 


This set of equations can be numerically integrated to determine the overall 
heat rejection rate and the average fluid and tube wall (radiator fin root) 
temperatures as functions of distance down the tube. 

For laminar flow in the tube the local Nusselt number can be used taking 
into account the variable viscosity effects. At present the best Nusselt 
number representation would be to use that developed by Chen (reference 3) 
for the radiant flux boundary condition in terms of the uniform wall flux 
Nusselt number, i.e., 

Nu(Gz)~ jo. 9ze Lp — 3 ) J Ni/f (Gz) ( 9 ) 

The variable viscosity effects should be well represented by using the 
results of Shannon and Depew (reference 5) to modify the constant property 
uniform wall flux Nusselt number Nu , i.e., 

q 

^ W (6z) 

Mu (Gz) = £ 0 92S + sin Lf ± Jj fJL J Nu^(Gz) no) 
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The viscosity ratio exponent as a function of Graetz number m(Gz) is given 
in reference 5 and the uniform wall heat flux Nusselt number as a function 
of Graetz number Nu^(Gz) is given in reference 2. 


Having calculated the tube wall and average fluid temperatures as a function 
of distance down the tube, the pressure drop may be calculated including 
the variable viscosity effects from 


dp_ __ __ /_£& W/J f£ I n ( (3z '> 

c/x TTfD+ (>*,/ 


where the viscosity ratio exponent as a function of the Graetz number 

n(Gz) is given in reference 5 and the viscosities are evaluated at the local 

temperatures (i.e. , at T(x) and T (x)) . 

R 


Due to the lack of data correlations for the local Nusselt number the 
calculations for turbulent flow are less refined that those for laminar flow. 
The overall heat rejection rate and tube wall and average fluid temperature 
variation with distance along the tube may be determined using the average 
Nusselt number given by equation (6) . The pressure drop may be calculated 
from 


it - 



( 12 ) 


where f (Re) is taken from the "Moody Diagram." 

An actual liquid loop radiator system involves much more than the single 
tube radiator considered in the preceding paragraphs. However the proper 
treatment of the fluid flow in a single tube is critical to the overall 
radiator system performance. Under nominal heat load conditions the heat 
transfer and pressure drop considerations for the radiator tubes can be of 
primary importance. In order to achieve maximum flexibility radiator 
designers strive to achieve as large of a heat load range as possible while 
maintaining temperature control. Several techniques are used for controlling 
the temperature at low heat load conditions, among them are: 
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o Bypass Control 

A value is used to bypass the fluid around the radiator at low heat 
load conditions. This type of control is usually unacceptable due to 
the possible freeze up of the radiator. 

o Regenerative Control 

Full flow is maintained through the radiative and a regenerative heat 
exchanger is used to control the fluid temperature. This technique 
generally results in very large pressure drops at low heat loads and 
is limited by the pump capabilities. 

o Selective Stagnation Radiators 

This technique uses parallel flow passages in the radiator panel and 
allows the fluid to stagnate (freeze) in selected passages as the 
heat load is reduced. The successful implementation of this 
technique depends on the accurate prediction of the pressure drop in 
each flow passage. Flow instabilities may arise which cause progressive 
freeze up of the panel. 

All of these control techniques depend on the accurate assessment of the 
pressure drop through the radiator tubes. This is a critical problem for 
fluids that are characterized by a pour point temperature. Since their 
viscosities are extremely temperature dependent, fluid temperature varia- 
tions in both the radial and axial directions have large nonlinear effects 
on the pressure drop. These effects are accentuated by the fact that these 
fluids are in the laminar flow regime at the low heat load conditions. 

Fluids that are characterized by a freezing point temperature usually have 
viscosities that are not strongly temperature dependent. Consequently, the 
variable viscosity effects are less critical for these fluids especially 
if they are in the turbulent flow regime. However, the use of these fluids 
in selective stagnation radiators requires the freezing and thawing 
phenomena to be understood in relationship to possible flow instabilities in 
the radiator panels. 
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Free convection effects are absent for liquid loop radiators operating 
under zero-g conditions. However, the design verification tests are 
conducted in a 1-g environment where free convection may significantly 
affect the radiator performance. 
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IV. 4 THERMAL SCALE MODELING CRITERIA 


The thermal scale modeling criteria for liquid loop space radiators may be 
developed from the energy equation, momentum equation and the boundary 
conditions. For steady state conditions and constant fluid properties the 
energy equation (3) may be written as: 

RePr V* - v/Q' <i3) 


the momentum equation (4) as 

■+ (14 > 


and the boundary conditions, equations (7) and (8), as 

-PePr-B- 4Nu (§'-&'*) 


(15) 


and 


4NU (S'- & } = ~ 3 (Q ’«' 


(16) 


Thermal similitude with geometric scaling then requires the following 
parameters to be preserved in the scale model: 

Re , Pr, 

If free convenction effects are present then the Grashof number Gr must also 
be preserved in the scale model. For the variable viscosity effects to 
be preserved in the scale model the (p/p) term (see equation 4) must be 
preserved. For a quasi real fluid (a fluid for which the viscosity is 
exponentially dependent on temperature, see reference 5) this term may 
be written as 


V 



(17) 
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Consequently, since most real fluid viscosities may be represented by this 
quasi real fluid over given temperature ranges, the preservation of the 
variable viscosity effects requires the preservation of the wall to bulk 
viscosity ratio. 


The thermal scale modeling criteria for liquid loop space radiators under 
steady state conditions can be summarized as 


£&*■= Pr* = 


'ecrtfpj-fjl 

4 > 


3 \ * 


= / 


(18) 


when free convection and variable viscosity effects are important the 
additional scaling criteria are 



(19) 


The detailed consideration of radiator freeze up for fluids characterized by 
a freezing point temperature is beyond the scope of the present study. 
However it is obvious that the nondimens ional freezing point temperature 
must be preserved in the scale model. Other factors involving transients 
may also be involved. 
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IV. 5.0 THERMAL SCALE MODELING TECHNIQUES 

The material preservation and temperature preservation thermal scale modeling 
techniques commonly used for radiation-conduction scaling are inappropriate 
for liquid loop radiators. The strong dependence of the fluid properties 
on temperature precludes the use of material preservation and it is extremely 
doubtful that a fluid could be found to exactly meet the scaling criteria 
for temperature preservation. Consequently, compromised scaling techniques 
must be considered. This is also true for scaling of systems involving 
gaseous convection (see Part II of this report) however, the scaling of 
liquid loop radiators is further complicated by the variable viscosity 
effects on heat transfer and pressure drop. 


The thermal scale modeling compromises to be used depend on the radiator 
characteristics of interest. Since flow instabilities and radiator flow 
stagnation are critical aspects for radiator concepts applicable to AAP 
and other advanced spacecraft, these aspects should be preserved as well as 
possible in scale model radiators. Various potential scaling techniques are 
developed in the following paragraphs . 

IV. 5 . 1 Fluid Preservation 

Since the fluid properties are critical to the radiator performance, it may 
be desirable to preserve the fluid in the scale model and use compromised 
scaling techniques to approximately preserve the temperatures. 

IV. 5. 1.1 Reynolds Number Preservation 

If the radiator operates with the coolant flow in different flow regimes 
for different heat load conditions, then the Reynolds number must be 
preserved in the scale model to preserve the operating conditions. Preser- 
vation of the Reynolds number may also be required for radiators with 

parallel flow paths in order to preserve the flow distribution pattern. 

% 
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IV. 5. 1.1.1 Radiator Fin Geometry Distortion 


If the radiator geometry is scaled geometrically except for the radiator fin 
width (tube spacing) which is kept full scale, then the radiator performance 
may be preserved in the scale model. 


The scale modeling requirements for this scaling technique are 

L*= I 

W* = Q* — (_* ( 20 ) 

AP*= (L*)~ z 

This technique preserves the fluid and radiator fin temperatures 

T*= ( X- Tp.)*- AT*= I (2D 


Since the temperatures are preserved the variable viscosity effects are also 
preserved 



( 22 ) 


however the free convection effects are not preserved since 




(24) 


This scaling technique has the advantage of preserving the radiator 
performance including the variable viscosity effects. Also the reduction of 
the Grashof number in the scale model may represent a way to simulate 
operation in a zero-g field. The disadvantage of this technique is the 
geometric distortion that doesn't allow reduction of the radiator width. 

IV.5.1.1,2 Radiator Tube Diameter Distortion 

If the radiator fin geometry is scaled geometrically approximate thermal 
similitude can be achieved in the model by a dis tori ton in the scaling 
of the tube diameter. Dietz and Fleming (reference 10) investigated the 
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use of this technique in preserving the fluid temperature change through 
the system. The scaling requirements for this case are: 

Re* = / 

p"= =• ( L *)* <25 > 

AP* •= (L*)~ S 

This technique preserves the radiator fin temperature and the fluid 
temperature change through the system 

Tz* = AT* = I (26) 


however the temperature difference between the fluid and tube wall is not 
preserved 


(T-r*f= L * 


(27) 


Consequently the variable viscosity effects are not preserved and the 
Grashof number is reduced 

Gr*— ( l*Y (28 

The major disadvantage of this technique is the large pressure drop require- 
ments for small scale sizes. 

IV. 5. 1.1. 3 Exact Geometric Scaling 

Exact geometric scaling results in reduced temperature changes in the fluid. 
The scaling requirements for this technique are 

/ 

W * = L* (29 

ap*-(l*)~ !! ' 
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This technique preserves the average radiator fin temperature 



(30) 


however the fluid temperatures are not preserved 

AT*= (T-Tg)*= L* (3D 

consequently the variable viscosity effects are not preserved and the Grashof 
number is reduced 

Gr*= (l*) 4 " ( 32 ) 

The similitude achieved with this technique is not quite as good as that for 
the radiator tube diameter distortion technique (IV . 5 „ 1 . 1 . 2) , however, the 
pressure drop requirements are much less severe. 


IV.5,1.2 Mass Flux Preservation 


If the Reynolds number is not preserved it is possible to preserve the fluid 
temperature change through the system while using exact geometric scaling. 

The scaling requirements for this technique are 

Re* - U* 

W*=r Q* - (U*) Z (33) 

' j 

^ p* =r ( if) ~ LAM1UAR. 

= f * TUZ0ULE_UT 

This technique preserves the average radiator fin temperature and the fluid 
temperature change through the system 

T z * — AT* = I (34) 
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however the temperature difference between the fluid and the tube wall is 
not preserved 


C r-T g f= 


J/* 

A/u* 


( 35 ) 


For fully developed laminar flow equation (35) can be written as 


(T-T*) 




(36) 


and for turbulent flow as 

(l* )* & 


(37) 


The temperature preservation is fairly good in turbulent flow, however, in 
laminar flow where variable viscosity effects can be important, the lack 
of temperature preservation precludes the preservation of variable viscosity 
effects. The major disadvantage of this scaling technique is the lack of 
Reynolds number preservation. 


IV. 5.1.3 Heat Transfer Coefficient Preservation 


Whereas the mass flux preservation scaling technique preserved the fluid 
temperature change through the system while using exact geometric scaling, 
the heat transfer coefficient preservation technique preserves the temperature 
difference between the fluid and the tube wall. 

IV. 5. 1.3.1 Turbulent Flow 


The heat transfer coefficient preservation scaling requirements for the 
turbulent flow case are 


iv* ■= 
<?* = 
A P* = 



(38) 
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This technique preserves the temperature difference between the fluid and tube 
wall 


Cr~T g f= I 09 ) 

However the fluid temperature change through the system is only approximately the 
same 

at *= (L*y°- zs 

The major disadvantage of using this technique for turbulent flow is that the 
lack of Reynolds number preservation may result in laminar flow for the scale 
model. 


IV.5.1.3.2 Laminar Flow 


The heat transfer coefficient preservation scaling technique is not applicable to 
fully developed laminar flow where the Nusselt number is independent of the 
Reynolds number. For developing laminar flow the average Nusselt number is 
generally taken as proportional to the cube root of the Graetz number. 




(41) 


For this case the scaling requirements are 


(L*f 

w* = (l*) 4 < 42 > 

q* - (L*y 

AP*= L* 

This technique preserves the temperature difference between the fluid and tube wall 

(~r-T*)*= I ( 43 ) 
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however the fluid temperature change through the system is not preserved 

AT* = (L*y* (44) 

The major disadvantage of using this technique for laminar flow is the non- 
preservation of fluid temperature change through the system. While this technique 
preserves the variable viscosity effects caused by radial temperature gradients 
those caused by axial temperature gradients are not preserved. 

IV. 5 . 2 Fluid Change 

IV. 5 . 2 .1 Temperature Preservation 

Use of the temperature preservation scaling technique depends on finding a 
suitable fluid for the scale model. Dietz and Fleming (reference 10) selected 
R-13B1 fluid for a 0.384 scale model of a space shuttle radiator panel which uses 
R-21 fluid. The slight difference in fluid Prandtl numbers neccessitated a 
slight distortion in the scaling of the tube diameter. Even though the freezing 
point of R-13B1 is -270^F as compared to -211%F for R-21, the results of their 
analytical investigation showed good thermal similitude over the entire range 
of radiator operation. The major disadvantage of this scaling technique is that 
the scale ratio is fixed at a given value depending on the model fluid. 

The use of this technique, in scale modeling a system in which variable viscosity 
effects are important, would probably be precluded by the lack of a suitable 
fluid for the model. 

IV. 5.2.2 Matched Viscosity Scaling Technique 

The thermal scale modeling techniques described in the preceding paragraphs are 
unsuited for laminar flow with variable viscosity effects with the exception of 
the radiator fin geometry distortion technique (IV. 5. 1.1.1) . This technique, 
however, does not allow the radiator fin width to be reduced in size. A solution 
to this problem might be found in a "matched viscosity scaling technique". This 
technique would use a fluid in the scale model that has a much larger viscosity 
than that of the prototype fluid at any given temperature and the prototype 
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performance would be simulated by operating the scale model at a higher tempera- 
ture. This would allow the radiator fin width to be scaled down in size. 


The scale model fluid would be chosen such that its viscosity dependence on 
nondimensional temperature matches that of the prototype f luid, i.e. , 

jj(* — ( POR O'* =■ I (4.5) 

This fluid selection fixes the characteristic model temperature since the model 
fluid viscosity at this temperature must equal that of the prototype fluid at 
the characteristic prototype temperature. 


The product of fin effectiveness and fin width is given by 

(ft L f ) * = ( T?) 3 d* ( 46 ) 

As can be seen from equation (46) the fluid thermal conductivity ratio has a 
strong effect on the scaling of the fin width and may counteract the reduction 
in size resulting from the increased model temperature. 


The local Nussett number, the nondimensional temperature difference between the 
fluid and tube wall and the nondimensional fluid temperature change through the 
system may be preserved by preserving the Graetz number. This can be 
accomplished by setting the model flow rate such that 


w 




(47) 


This technique geometrically scales the radiator except for the fin width which 
is distorted according to equation (46). Further study of this technique is 
required before its applicability can be ascertained. 
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IV. 6 CONCLUSIONS AND RECOMMENDATIONS 


Thermal scale modeling of liquid loop space radiators operating in the turbulent 
flow regime appears quite promising. Three scaling techniques may be used in 
one geometrically scaled model radiator configuration. These techniques use the 
same fluid as the prototype and preserve the radiator fin temperature. The heat 
rejection is scaled as the square of the scale ratio. The fluid temperature and 
temperature gradients vary somewhat between the scaling techniques. These 
techniques are implemented by adjusting the model flow rate as follows: 


Reynolds number preservation 


w* = L* 


Mass flux preservation 


w* = (L*V 


Heat transfer coefficient preservation w * = 


2.25 


The freeze up phenomena scaling requires further study. 


Thermal scale modeling of liquid loop space radiators operating in the laminar 
flow regime where variable viscosity effects are important is not as promising. 
One scaling technique that preserves the radiator performance does not allow a 
reduction in the radiator width dimension. Another technique, "Matched Viscosity 
Scaling Technique," shows some promise, however, it depends on the selection of a 
proper fluid for the scale model and further investigation of this technique is 
required . 

It is recommended that the thermal scale modeling study of liquid loop radiators 
be continued with emphasis on the high load ratio radiator systems required for 
AAP and other advanced spacecraft. The results of this study should then be used 
to experimentally verify and demonstrate the use of thermal scale modeling in the 
design, development and test verification of liquid loop radiator systems. 
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PART V 

EVALUATION OF EXISTING NASA/MSC FACILITIES 
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V.l INTRODUCTION 

Thermal verification tests of Apollo Application Program spacecraft will be both 
demanding and difficult • The physical size and complexity of the AAP configurations 
are such that full scale testing of the assembled vehicles is not possible utilizing 
existing thermal /vacuum test facilities. 

Consequently, scaled thermal model testing as an alternative to full scale testing 
must be considered. Scaled model tests may provide certain advantages in terms of 
reducing test facility and test model costs. 

V. 2 AAP CONFIGURATIONS 

The Apollo Applications Program is primarily concerned with the Skylab Spacecraft 
configuration as shown in Figure 1. The figure shows the completely assembled 
spacecraft with a docked command and service module. The physical complexity of 
the assembled spacecraft and the potential for significant solar shadowing are 
apparent from the figure. 

In addition to the AAP-Skylab other large manned spacecraft presently being studied 
include the Space Station, the Space Base and the Space Shuttle. 

All of these spacecraft configurations present similar thermal simulation test 
problems because of their size and physical arrangement. Practical thermal 
testing of the assembled spacecraft can probably only be accomplished with scaled 
models . 

V. 3 EXISTING NASA/MSC FACILITIES 

Existing NASA/MSC facilities for thermal /vacuum testing are located in both the 
Space Environment Simulation Laboratory and the Space Environment Effects Labora- 
tory. The Space Environment Simulation Laboratory consists of chambers A & B. The 
Space Environment Effects Laboratory consists of chambers D & E and six smaller 
chambers. A description of the capabilities and characteristics of all of the 
space environment test chambers is included in reference 1. 
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V. 4 FACILITY REQUIREMENTS FOR SCALE MODEL THERMAL TESTS 

Facilities required for spacecraft thermal performance tests consist of vacuum 
chambers with a cold wall to simulate deep space and some form of thermal simulation. 
Solar thermal energy simulation requires simulation in terms of total energy, spec- 
tral distribution and collimation. Planetary thermal emission and albedo simulation 
are usually not required except for special tests. 

General requirements of test facilities for scaled model thermal tests relate to 
vacuum chamber size, pressure and thermal environment. The required chamber size 
is determined by the scale model size. Recent research (reference 2) has indicated 
that scaled models smaller than 1/10 prototype size incur unacceptable temperature 
errors, while scaled models smaller than 1/ 6-1/8 prototype size are increasingly 
error sensitive to small changes in model scale ratio. This indicates that thermal 
testing of the Sky lab cluster would require a scaled thermal model approximately 
12-20 feet in length and 4-6 feet in diameter, not including solar arrays. While 
the Space Shuttle Orbiter would require a scaled model 15-20 feet in length and 
7-9 feet wide. 

Vacuum chamber pressure requirements for scaled thermal models will not differ 
from the requirements for full scale thermal model tests. 

Proper thermal environment simulation is an important requirement for scaled thermal 
model tests. Both the solar simulator and the deep space simulator cold wall must 
effectively simulate the actual thermal environment. The actual cold wall temper- 
ature and emittance are not any more critical than what is required for full size 
thermal models. However, large temperature gradients in the cold wall or large 
variations in the wall emittance could have a significant though local effect for 
a case where the model and chamber are approximately the same size. 

Solar simulator errors in intensity, uniformity, collimation and spectral match 
can all seriously affect the results of thermal tests. However, these errors do 
not affect scaled models anymore than full size models. However, scaled models 
tested in large chambers will normally utilize only a small section of the solar 
simulator which may have different characteristics than the overall simulator. 
Consequently, thermal simulator characteristics must be accurately determined 
for all thermal tests and particularly so for scaled model tests. 
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V.5 EVALUATION OF NASA/MSC FACILITIES IN COMPARISON TO REQUIREMENTS 
There are four existing NASA/MSC facilities which can be considered for use in 
thermal scaled model testing. Chambers A & B in the Space Environment Simulation 
Laboratory and chambers D & E in the Space Environmental Effects Laboratory, The 
remaining chambers in the Space Environmental Effects Laboratory (FGHIK & N) are 
either limited in size, used for special functions or do not have adequate thermal 
simulation. 

Table I shows the interior working dimensions and solar simulator dimensions for 
chambers A, B, D & E. 



A 

TABLE I 
CHAMBER 
B 

D 

E 

Working 

75 ’L 

27 'L 

15 ’L 

10 'L 

Dimensions 

25 ’D 

13’D 

6'D 

4'D 

Solar 

Top 13' D 

7.5'D 

3 . 5 'D 

3'D 

Simulator 

Dimension 

Side 13 'x33 ' 





The approximate dimensions of the Sky lab cluster are 120 1 x 40 V. Space Shuttle 
orbiter dimensions are approximately 120 1 x 60 1 while the assembled orbiter/ 
booster is approximately 200 1 x 140*. Table II shows the various size scaled 
models of Skylab and Space Shuttle that could be tested in the NASA/MSC chambers. 

This table indicates that chamber "A" can be used to test any realistic size 
scaled thermal model of either Skylab or Space Shuttle. Chamber B however, can 
be used for testing these models only with modifications to provide more extensive 
thermal simulation capability. Chamber D with modifications to the thermal simu- 
lation capability, can be used for tests of small models of Skylab or Space Shuttle. 
Chamber E can be used only for tests of very small sized models, for example, a 
1/12 scale model of Skylab or the Space Shuttle orbiter would be necessary. Such 
small models would incur probable temperature errors of approximately 6-8 percent 
of absolute temperature. Chamber E could however, be used to test a small 
scaled model of the Skylab cluster minus the CSM and ATM. 
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TABLE II 
CHAMBER 


Conf ig. 

A 

B 

D 

E 

SKYLAB 





1/4 

X 

(X) 



1/6 

X 

(X) 



1/8 

X 

(X) 

(X) 


1/10 

X 

X 

(X) 

* 

SPACE SHUTTLE 

ORBITER 




1/4 

(X) 




1/6 

X 

(X) 



1/8 

X 

(X) 



1/10 

X 

(X) 

(X) 



X - Test with no mod 

(X) - Test with mod to chamber 

* - Skylab Test minus CSM & ATM 


V. 6 CONCLUSIONS 

Thermal scale model tests of reasonable sized models (1/4-1/10) of both Skylab 
and Space Shuttle Orbiter can readily be performed in NASA/MS C chambers A & B, 
Chamber A modifications to the thermal simulator are required only in the case 
of a large (1/4) model of the Space Shuttle Orbiter. Extension of the thermal 
simulator capability of Chamber B will be required for most proposed model tests. 

A significant saving in test costs can probably be made by using chamber D where 
possible. Table II indicates that only quite small models of either Skylab or 
Space Shuttle can be tested in this facility and then only with major modifications 
to the thermal simulation capability. However it is possible that wider use could 
be made of chamber D by confining thermal model testing to include only the areas 
of the vehicle which have critical temperature control requirements. 
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